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Evolution of a Computer Application 

John J. Wavrik 
Department of Mathematics 
University of California - San Diego 


Abstract 

In this article I discuss the process of producing a computer software 
system for mathematical research or instruction. I show how a 
mathematician can create a special-purpose computer language to fa¬ 
cilitate development of a system. 

Mathematical work benefits from the use of software that can be modi¬ 
fied and extended as it is used. I show the process of creating a soft¬ 
ware system of this type by discussing several critical stages in the 
development of Groups32, a system for working with groups of low or¬ 
der. 
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A version of the Groups32 system being discussed in this article is accessible on 
the Internet at http://math.ucsd.edu/~iwavrik/q32/ . This web page contains in¬ 
structions for accessing Groups32 via telnet and some written material including 
a sample session. 
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Introduction 


I address this article to mathematicians 
and mathematics educators. It is not about 
the use of a piece of software but rather 
about a methodology for producing soft¬ 
ware and about the electronic representa¬ 
tion and manipulation of mathematics. In 
the article I show how a mathematician 
can produce a special purpose language 
for a subject area to facilitate writing pro¬ 
grams in that area. 


Computer programs are often "end prod¬ 
ucts" -- like books, they remain relatively 
fixed after creation. The design and im¬ 
plementation of large software projects are 
most often in the hands of teams of ex¬ 
perts in computer programming. The end 
user is not interested in the programming 
language or the developmental techniques, 
and, in fact, these things become invisible 
in the end product. The choice of lan¬ 
guage, however, can be a significant factor 
for ease in development, maintenance, 
and modification. 
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In this article I discuss production of dy¬ 
namic software systems that are targeted 
at specialized research and instructional 
areas and are extended and modified as 
part of their use. The user, designer, and 
implementer are often the same person. 
We need, for this work, a language and 
methodology that allow software to be de¬ 
veloped and modified easily. We need to 
make it feasible for a mathematician, 
whose expertise lies in other areas, to be¬ 
come involved in the task of creating lan¬ 
guages and software systems. 


A source of difficulty in programming with 
conventional general-purpose languages is 
the gap between the means of expression 
provided by the language and the concepts 
of the subject area. We can reduce this 
gap by using a language specifically de¬ 
signed for a particular subject area. To do 
this in the realm of conventional compiled 
languages would require writing a custom 
compiler - a formidable task. I will discuss 
and illustrate an approach that is the 
moral equivalent of writing a compiler - 
but that is viable for a working mathemati¬ 
cian. 
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It has been clear for some time that one 
cannot portray the process of creation by 
writing an article that describes it. The 
reader must have direct experience with 
the systems under discussion. This journal 
facilitates exposition that combines written 
material and access to software. We are 
interested in electronic representation and 
manipulation of mathematics and in fac¬ 
tors that influence the development of 
software. The best way to address these 
matters is to discuss the development and 
evolution of a particular software system. 
I have chosen a system written for ele¬ 
mentary group theory, Groups32, because 
the subject matter is familiar to most 
mathematicians. 


I will not distinguish in this article between 
research and instructional systems. 
Groups32 is not an "instructional program" 
- that is, it does not contain expository 
material to teach group theory. Rather, it 
is a research program that has been made 
suitable for instructional use. It provides 
students with a tool for experimentation. 
In a theoretical abstract algebra course, I 
have found it a good system for providing 
experience, examples, and experimenta¬ 
tion. 
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History of Groups32 

My first step in developing what came to be Groups32 was inspired by 
Kenneth Almquist's electronic publication in 1989 of a set of tables for 
groups of orders 1-16. I wrote the first stages of the program (then 
Groupsl6) just to solve a problem stemming from an accident: Alm¬ 
quist's posting contained some tables that represented isomorphic 
groups. A natural question was to determine what one must know 
about groups of order 1-16 to determine quickly whether two tables 
represent isomorphic groups. Eventually, Groups32 evolved into a 
more substantial program that now contains tables for all groups of 
orders 1-32 and a collection of algorithms to compute a variety of in¬ 
formation about these groups. I developed the program using the 
same techniques as for several special-purpose research systems. 


Table of Contents 
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Tips on Reading This Article 


This article is intended to be read in con¬ 
junction with a computer application that 
is accessible on the Internet. 

1. The Grouos32 Web page contains in¬ 
structions for accessing Groups32 via 
Telnet and also includes some written 
material. 

I use this site with students. This way of 
making software available has the advan¬ 
tage that students can use it at home with 
any computer or operating system. The in¬ 
structor can add material during the 
course. The Telnet version is not pro¬ 
grammable. It does not require knowledge 
of a programming language. 


2. I include with this article, for readers 
who would like to experiment with 
programming, a ZIP file containing ex¬ 
ecutable versions of an early (1990) 
version of the groups program and the 
current version. These executables will 
run only under Windows 95 or later. 

a. Groupsl6 was written in 1990. I include 
it for paleontological purposes. 

b. Groups32 is the current version. I have 
used it in courses in abstract algebra at 
the undergraduate and beginning gradu¬ 
ate level and for REU projects. 
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Unless you are unusually well endowed - 
with extremely good eyesight and a very 
large monitor — it is neither realistic nor 
necessary to have both the article and the 
software visible at the same time. It is 
enough to have them running in separate 
windows that you can switch between 
quickly. 

A. Start the Groups32 software (either 
via Telnet or using a downloaded ver¬ 
sion on your own computer). Run the 
software in a large window. 

B. Start the article (using Adobe Reader) 
in a separate large window. The article 
is intended to be read without scroll¬ 
ing. Its pages fit comfortably within a 


large Adobe Reader window. Make 
sure the view is set to "fit to width" 
and the window is as large as possible. 
Use the navigation arrows on the tool¬ 
bar at top and navigation links in the 
article to move from screen to screen. 
If you are reading with the Adobe 
Reader inside a browser window, you 
may need to hide some of the browser 
toolbars. 

C. Switch back and forth by using what¬ 
ever method your operating system 
provides for switching between win¬ 
dows. (In the case of Windows you 
can switch using the buttons on the 
taskbar at the bottom of the screen.) 
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FORTH : 

I have used the computer language Forth as the base language for this 
work because it is well suited for the task of implementing special- 
purpose languages and creating highly interactive mathematical sys¬ 
tems. 

It is not the purpose of this article to make readers expert in Forth. 
The interested reader will find a variety of resources (Forth implemen¬ 
tations, books, tutorials, etc.). Most of this article requires the devel¬ 
opment of a reading knowledge of the language. This means that you 
should be able to understand the simpler code and get the flavor of 
the more technical code. It also means that, if you are using one of 
the downloadable versions, you should be able to write some simple 
procedures. You should be able to read lightly through the main parts 
of the article first and then look at the section on Forth to get back¬ 
ground for more detailed reading. 
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To program in Forth you need a development system. The downloadable versions of the Groups 
system include Win32Forth, a high quality Forth development system written by Tom Zimmer and 
Andrew McKewan, which is in the public domain. It is distributed with permission of the authors. 

I would like to express my thanks to Tom Zimmer for making it possible to provide a develop¬ 
ment system with this article. 


Table of Contents 
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For the benefit of those who are trying to execute some of the exam¬ 
ples in the articles on the downloadable systems, we have color-coded 
the code fragments that appear in the article. 

The Color Coding Explained 

• red code is found in the 1990 version, Groupsl6. 

• green code is found in the current Groups32. 

• blue represents code that was not in the 1990 version and 
did not survive to the current Groups32. For the most part, 
this code represents important stages in the transition. 
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The next section is a brief introduction to Forth. If you are using one of 
the programmable versions, you should try some of the examples in 
this section and conduct experiments to get the "hang" of the lan¬ 
guage. 

The section after that discusses an early version of the groups pro¬ 
gram. The most important documentation in Forth is a glossary. 
Appendix A (a separate PDF file) is a glossary of the standard Forth 
words used in this article, and Appendix B is a glossary of the words de¬ 
fined in Group32. Individual words in the glossaries are bookmarked. A 
description of an individual word can be located most rapidly if the 
Adobe Reader is set to display bookmarks when using a glossary. 

I suggest that you run the Groupsl6 (1990) program in one window 
while reading the article in another, and try some of the examples. 
The glossary should also be available for reference - it can either be 
open in another window or printed. 

Table of Contents 
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The Forth Language 


Forth is a language designed to create lan¬ 
guages. It is easy to learn by understand¬ 
ing how it works! Unlike most conventional 
languages, its syntax is based on seman¬ 
tics: The emphasis is on what is done 
rather than on how it is said. The quickest 
way to introduce you to the process of 
programming is to examine and discuss 
some code, and that is what I will do in 
the remainder of the article. I start with 
some background information. 


The Dictionary 

The Forth language is a collection of com¬ 
mands, which we may think of as words in 
a dictionary. The name for a command 
(which is called a "word") can be any col¬ 
lection of printable characters delimited by 
spaces. Each command has an associated 
action. The dictionary consists of the 
names and some stored representation of 
the action. 
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The process of programming in Forth con¬ 
sists of adding new words (and actions) to 
the dictionary. A glossary of standard 
Forth words used in this article is found in 
Appendix A . Individual words can be ac¬ 
cessed using the bookmarks in Adobe 
Reader. 


Forth maintains its dictionary in chrono¬ 
logical (rather than alphabetical) order. 
The dictionary looks like a linked list: 
When a new word is defined, it is added at 
the end of the dictionary. 


Previous Name 


Action 


Name 


Action 
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Interpreting 

The Forth system has two modes. The first 
- interpreting - occurs when the system 
reads a line of input. When a word is en¬ 
countered in the input stream, the diction¬ 
ary is searched, starting from the most re¬ 
cent entry. If the word is found, its action 
is performed. If a word is not found, the 
system tries to interpret it as a number. 

Since a word's action is performed when 
the word is encountered, any parameters 
it needs for its action must be already 


available. Thus the words (or numbers) 
that are used as parameters must occur 
before the word in the input stream. Some 
means must be provided for storing these 
parameters - a stack is used for this. By 
convention, most words act by removing 
their arguments from the stack and placing 
their results on the stack. (While numbers 
are not stored in the dictionary, they can 
be thought of as words whose action is to 
put themselves on the stack.) 
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Example: 

3 2 + 

This sequence (reading left to right) puts 3 
on the stack, and then puts 2 on the stack. 
+ is a word that takes two numbers from 
the stack, adds them, and puts the result 
(5 in this case) on the stack. Here are 
some examples with two operations: 


3 2 + 7* leaves 35 on the stack 
3 2 7 *+ leaves 17 on the stack 
3 2 7 + * leaves 27 on the stack 

There is no hierarchy of operations or use 
of parentheses to postpone operations. 
The sequence 3 2 7 puts these numbers on 
the stack (with 7 on top). The operations 
+ and * are performed on the top two 
numbers on the stack. 


NOTE: I f you try these examples you will not see the results. The numbers on the 
stack are not usually displayed. Read on to find out how to see numbers on the 
stack. 
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The stack is a central feature of Forth. It 
(rather than named variables) is the 
means by which words communicate with 
one another. In most languages data items 
must be named before they can be used. A 
typical Forth word takes some unnamed 
parameters from the stack, performs some 
operation on them, and leaves the re¬ 
sults) on the stack. The action of a word 
is usually documented by a "stack dia¬ 
gram" that shows what the word expects 
on the stack before it acts and what it 


leaves after. Often the programmer adds a 
further note of explanation, but sometimes 
this is obvious from the name chosen for 
the word. The stack diagram for + or * is 

(nl n2 - n3) 

which indicates that each expects two in¬ 
tegers on the stack and leaves one inte¬ 
ger. The names + and * indicate quite 
clearly what these words are intended to 
do. 
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Seeing what is on the stack 

If you type 3 2 + (and press enter) you operations that you must know about for 
will not see anything! The result is on the seeing results: 
stack, which is not visible. There are two 


• 

(n- ) 

This is a dot (period). It prints the num¬ 
ber on top of the stack (and removes it). 

.s 

(--) 

This (dot S) prints the stack without re¬ 
moving anything. The numbers are 
printed horizontally with the top of the 
stack on the right. 
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Example: 

3 2 7 + .S [2] 3 9 .. 

We started with 3, 2 and 7 on the stack, but the 
top two numbers were added leaving 3 and 9. In 
Win32Forth the number of elements on the stack 
is shown (by .S) in square brackets followed by 
the stack elements themselves. Win32Forth also 
prints dots at the end of each line showing how 
many numbers are on the stack. In this case two 
dots indicate two numbers on the stack. 

* .S [1] 27 . 

After multiplication, .S shows that there is only 1 
number on the stack and it is 27. In both cases 
.S does not remove numbers. 

. 27 

Finally we type . (dot) and the 27 is printed and 
removed from the stack. 
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There are several words in the Forth language whose action is to rearrange the items on 
the stack. Here are the basic ones: 


DUP 

( n - n n ) 

duplicate number on top of 
stack 

DROP 

(n-) 

drop number on top of stack 

SWAP 

( nl n2 - n2 nl ) 

swap numbers on top of stack 

ROT 

( nl n2 n3 - n2 n3 nl ) 

bring 3 rd number to top 

OVER 

( nl n2 - nl n2 nl ) 

copy second element to top 


To see the number on top of the stack without removing it, you can use dup. (DUP fol¬ 
lowed by dot). This will duplicate the number before printing it. 
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Compiling 

Programming in Forth is equivalent to add¬ 
ing new words to the dictionary. Here is 
the main way new words are added. The 
Forth word : (colon) is used to start a 
definition - it is followed by the name of 
the word to be added. The definition is 
terminated by the Forth word ; (semico¬ 
lon). The action of the : (colon) is to cre¬ 
ate a new dictionary header and switch the 
system to compiling mode. The action of ; 
(semicolon) is to compile a termination 
word and to return the system to the in¬ 
terpreting mode. The other words in the 
definition describe a sequence of actions to 
be taken when the new word executes. At 
the time of definition, these words do not 


carry out their action - the action is stored 
in the new definition. 1 

Example: 

: 2* DUP + ; 

We have chosen to name the new word 
2*. The action of this word is to perform 
the action of DUP followed by the action of 
+ . This adds the number on the stack to 
itself. 


1 For the curious: Traditional Forth systems store the 
address (within the dictionary) of the component words. 

In some recent versions of Forth, words are subroutines 
in machine language and component words are compiled 
as subroutine calls. In either case, no dictionary search 
for component words is required when a word executes - 
the dictionary search takes place during compilation. 
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Compiling the action of component words 
results in a major speed increase without 
sacrificing the interactive nature of the 
language. The most time-consuming step 
is the search of the dictionary. The new 
word AddTest: 

: AddTest 32+ DROP ; 

executes 10,000,000 times in 1 second. 
The sequence of component words: 

32+ DROP 

takes 600 seconds to evaluate 10,000,000 
times. This increase of 600x in execution 
speed occurs because, in interpreting 
mode, the dictionary must be searched 


each time the sequence is executed. In the 
compiled word, the dictionary is searched 
only when AddTest is compiled. Compila¬ 
tion is incremental: It occurs only when a 
new word is added and affects only the 
new definition. (It is neither desirable nor 
necessary to recompile an entire system to 
add a new command. 2 ) Indeed, new defi¬ 
nitions can be made and compiled while in 
the midst of an interactive session. Execu¬ 
tion speed for (compiled) Forth words is 
comparable to the execution speed for 
conventional compiled languages. 


2 Modern computer algebra systems (such as Maple 
and Mathematica) use a similar approach: They provide 
an extensible, interactive environment while retaining 
execution speed by compiling new procedures to an in¬ 
ternal code. 
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Terminology and Naming Convention 

The words : and ; function as if they were called something like "DEFINE" and "END-DEF". 3 
The definition of 2* might seem clearer if we were to write 


DEFINE 2* 

Make a new dictionary header with name 2* 

DUP + 

body - describe the action of the new word 

END-DEF 

end the definition 


3 


One can, in fact, easily add these new names to the dictionary and use them instead. 
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Forth tends to use punctuation marks as 
names for frequently performed opera¬ 
tions. Typing speed is not the only rea¬ 
son: These punctuation marks are often 
used as part of a naming convention to in¬ 
dicate what new words do. The words : 
and ; are often used by programmers to 
suggest opening and closing something: 
Win32Forth has words COMMENT: and 
COMMENT; that are used to open and 
close multi-line comments. 

The word . (dot), used to print the number 
on the top of the stack, is often used as 
part of the name of a word that prints 


something. A Forth programmer who sees 
a word with a dot in it will automatically 
assume that this word is intended to print. 
Later in this paper, for example, we will 
define a word "Set." that (as the naming 
convention suggests) is a word that prints 
a set. 

It is extremely helpful, when writing Forth 
programs, to choose names well. A good 
name should suggest its action. The use 
of programming conventions involving the 
punctuation marks is also very helpful. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 26 


Control Structures 

Forth provides control structures, such as 
IF .. THEN, which are like those found in 
most languages. IF and THEN are Forth 
words that are in the dictionary and, like 
other Forth words, they have associated 
actions. As noted already, most Forth 
words do not act during compilation - they 
are, instead, compiled into the current 
definition. Control structures change the 
flow of execution and are implemented by 
words that carry out their action during 
compilation. The dictionary header for a 
word contains both the name of the word 
and some extra information about it. A 
word can be tagged (a so-called "immedi¬ 
ate" word) so that its action is performed 
immediately, regardless of the state of the 


system. Because these words make 
changes in a dictionary entry, they are 
used only during compilation. 
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The placement of words such as IF and 
THEN in sentences can be best understood 
if one knows their action. Let A, B and C 
be collections of words. We stipulate that 
A must leave a number on the stack - zero 
is interpreted as FALSE, any non-zero 
number as TRUE. We do not want to exe¬ 
cute A, B, C, in sequence - we wish to 
execute B only if the result of A is TRUE. A 
flow diagram for this behavior is: 


Before B there must be a test of the num¬ 
ber on the stack. If this number is 0 
(FALSE), then we want to jump past B to 
the point right before C. Our source code 
must indicate where the test takes place 
and the target of the jump. The word IF 
compiles the conditional branching instruc¬ 
tion, and the word THEN indicates the tar¬ 
get. 
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Thus, in the definition of a new word, we This will be compiled into the dictionary 
write: entry of the word to look like this: 

A IF B THEN C 


code for A 

7BRANCH 

target_addr 

code for B 

code for C 


The target address is right before the code 
for C. IF compiles the conditional branch¬ 
ing instruction and leaves space for the 
target address (which it, of course, does 
not know at the moment), so it leaves on 
the stack the location of the missing ad¬ 
dress. The word THEN occurs after B and 
THEN knows where it is! It also knows 
(because IF left the information) where it 
must fill in the address of where it is. So 
THEN fills in the proper target address - 


and we are left with compiled code that 
looks like the diagram above. 

Notice that IF and THEN do not themselves 
occur in the dictionary entry of the word 
currently being defined (although we 
sometimes talk as if they do). Their job is 
to act during compilation. Their action is to 
leave behind code that will, when it exe¬ 
cutes, provides the proper control flow. 
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The conditional branching instruction fol¬ 
lows the convention of other Forth words: 
It removes the number tested (its argu¬ 
ment) from the stack. 

Example: 

: ABS 

DUP 0< IF NEGATE THEN ; 

This word has the stack diagram ( n - n' ) 
where n' is the absolute value of n. The 
word 0< has the stack diagram ( n - t/f ) 
and NEGATE has the diagram ( n - n' ) 
with n' = -n. 


ABS negates a negative number but not a 
non-negative number. The definition can 
be written up with the stack comment and 
additional comments: 

: ABS ( n - n' ) 

\ n' is the negative of n 
DUP 0< IF NEGATE THEN ; 
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Besides IF .. THEN, there are 
structures: 

IF..ELSE..THEN 
BEGIN..UNTIL 
BEGIN..WHILE..REPEAT 
DO..LOOP 


other control Forth does not have a separate program 

called "the compiler" - the functionality of 
a conventional compiler is delegated to the 
action of the "immediate" words. 
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Data 

The last thing that must be discussed in 
this brief introduction is how data is han¬ 
dled. The word : (colon) is called a defin¬ 
ing word because it is used to add new 
words to the dictionary. There are other 
defining words that come with a Forth sys¬ 
tem, and Forth provides a mechanism for 
the user to make others. The defining 
words for data items create new words 
(children) in the dictionary and allocate 
memory for the data associated with those 
words. They also specify the run-time ac¬ 
tion of their children. The child word usu¬ 
ally will put either the address of the data 
item or the item itself on the stack. For in¬ 
teger data, we use three defining words: 
CONSTANT, VALUE, and VARIABLE. Here is 


how each of these creates a new word in 
the dictionary: 
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Example 

When word is created 

When word is used 

5 CONSTANT xxx 

Create a constant xxx 
initialized to 5 

xxx puts 5 on the stack. A 
constant cannot be 
changed once it is de¬ 
fined. 

5 VALUE yyy 

Create a value yyy ini¬ 
tialized to 5 

yyy puts 5 on the stack 
but 7 TO yyy will change 
the stored value to 7. 

VARIABLE zzz 

Create a variable zzz 
(some implementa¬ 
tions initialize to 0) 

zzz puts the address of 
the storage on the stack. 
Words @ (fetch) and ! 
(store) can access and 
store data. 


Table of Contents 
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Groups'! 6 - 1990 version 


The numbering of groups in this section is 
different from that in later sections and in 
the current Groups32 program. (Some of 
the groups in this section are isomorphic.) 


Groupsl6 is an early version of the groups 
system. The system originated with a 
posting to Usenet in 1989 by Kenneth 
Almquist. He had generated tables for 
groups of orders 1-16, apparently as an 
exercise in combinatorics programming. In 
his original posting some of the tables rep¬ 
resented isomorphic groups. Groupsl6 
started in an effort to find out how much 
(or how little) one must know about 


groups to distinguish which tables were 
isomorphic. 

Data Representation 

Input and Output 
Some First Procedures 

Evolution of ORDERS 

Finding Isomorphisms 

Return to Table of Contents 
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The Almquist tables looked like this: 

1 group of order 1: 

A| A 

1 group of order 2: 


A| A B 
B| B A 


1 group of order 3: 


A| A B C 
B| B C A 
C| C A B 


2 groups of order 4: 


A| A 

B 

C 

D 

A| A 

B 

C 

D 

B| B 

A 

D 

C 

B| B 

C 

D 

A 

C| C 

D 

A 

B 

C| C 

D 

A 

B 

D| D 

C 

B 

A 

D| D 

A 

B 

C 


1 group of order 5: 


A| A B C D E 
B| B C D E A 
C| C D E A B 
D| D E A B C 
E| E A B C D 
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An interesting thing happened: Some 
readers pointed out that Almquist had 
posted more groups of certain orders than 
the literature claimed. He traced the prob¬ 
lem to the fact that his algorithm had not 
detected that some of the original tables 
were isomorphic. A subsequent posting 
corrected the error. However, this raised 
an interesting question: How much do you 
need to know about a group of order 1-16 
to be able to easily detect isomorphism? 


It is obvious that, even for groups of order 
16, one does not try all possible bijective 
mappings to see if two groups are isomor¬ 
phic. Instead, one hopes to look at prop¬ 
erties of the groups that must be pre¬ 
served under isomorphism. The original 
steps in the system were to represent the 
data in a suitable way, and to implement 
words that computed various information 
about the groups. 
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Data Representation 


The original version was implemented on a 
computer that had only 64K of memory. 
The version of Forth used 16-bit integers 
and so could address only this much mem¬ 
ory. The Forth system and DOS used 26K 
of memory - and a nicer editor had been 
added to bring this to 44K. This left just 
20K for the group theory system and data. 

Fortunately the Forth systems commonly 
available at that time provided a rather 
clever way to use the hard disk for virtual 
memory. The disk was organized into 
1024 byte blocks, and the system provided 
a word BLOCK with stack diagram ( n - 
addr ). Using 5 BLOCK, for example, the 
system would check if block 5 on the disk 


had already been loaded into memory. If 
not, it would be read from disk into a 
buffer. In either case, the base address of 
the buffer would be returned. With this 
scheme, one could create the illusion that 
a large amount of data was in memory. 

For purposes of this article, I have moved 
the earlier code to a modern (32-bit) ver¬ 
sion of Forth on a computer that has 128 
megabytes of memory. I therefore put the 
tables into a section of the dictionary. To 
make the earlier code run without rewrit¬ 
ing, I redefined the word BLOCK to refer to 
parts of this section of memory. 
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The original tables were converted to an internal format by stripping everything except the 
interior squares of the tables. These were arranged sequentially by row: 

AABBAABCBCACABABCDBADCCDABDCBAABCDBCDACDABDABCABCDEBCDEACDEABDEA 

BCEABCDABCDEFBCAEFDCABFDEDFEACBEDFBACFEDCBAABCDEFBCDEFACDEFABDEF 

ABCEFABCDFABCDEABCDEFGBCDEFGACDEFGABDEFGABCEFGABCDFGABCDEGABCDEF 

ABCDEFGHBADCFEHGCDABGHEFDCBAHGFEEFGHABCDFEHGBADCGHEFCDABHGFEDCBA 


etc. 


This amalgamation of the group tables was 
put in Blocks 1-10. They were spaced so 
that tables did not straddle disk blocks. In 
order to locate a particular group, it is 
necessary to know the order of the group, 
the block in which the data are contained, 
and the number of bytes from the begin¬ 
ning of the block where the data are lo¬ 
cated. With this information we can calcu¬ 


late the address in memory for a particular 
table. 


Block 0 was used to store the triples 
Gnum, Gblk, Goffs for each group. 


0 CONSTANT IDXBLK 65 CONSTANT ID 


0 VALUE GNUM 0 VALUE GORD 

0 VALUE GBLK 0 VALUE GOFFS 
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At this point I made an important decision. 
Most of my computations would be for a 
certain group, fixed for the moment. 
Rather than write procedures to which I 
would have to pass the group number as a 
parameter, I designated the group of in¬ 
terest at the moment (given by number) 
as "the current group", and I wrote all pro¬ 
cedures to refer to that current group. 

My word to assign the current group is 
>group - it stores the group number as 
the value Gnum, as well as the triple men¬ 
tioned above. It is therefore not necessary 
to go back to block 0 each time we do a 
group computation, nor is it necessary to 
re-read the table from disk for each com¬ 


putation. Once the current group is speci¬ 
fied, vital information about it is read into 
memory. 

In the current version of Groups32, all the 
tables are stored in memory and the word 
>GROUP is replaced by an appropriate 
version. The virtual memory scheme is 
not used, and it is not essential for the 
reader to understand how it is used. It is, 
however, interesting that (1) there was 
such a mechanism for dealing with the lim¬ 
ited memory of early computers and (2) 
code written to use this scheme can be 
updated fairly easily. Click for more in¬ 
formation about virtual memory and the 
code for >GROUP. 
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: >GROUP ( n — ) \ set GORD GBLK and GOFFS 

DUP TO GNUM 

6 * IDXBLK BLOCK + DUP @ TO GORD 
2+ DUP @ TO GBLK 
2+ @ TO GOFFS ; 


[To run this on a 32-bit system using the original data file, the @ needs to be replaced by W@, which fetches a 16-bit value.] 


We now have the group data stored, to¬ 
gether with information needed to locate 
the start of a particular table. Next we 
want to implement multiplication - which 
means to locate an element in a particular 
row and column of the table. If we number 
the rows and columns from 0 to Gord-1, 
the entry in the i-th row and j-th column 
will be located i*Gord + j bytes from the 
start of the table. 


Remember that we are storing group ele¬ 
ments as letters A, B, C, ... that a com¬ 
puter is representing by ASCII codes 65, 
66, 67, ... . If we want to multiply B and C 
in the current group, we must first sub¬ 
tract 65 (the ASCII code for "A") to get the 
row and column, and then compute the 
location of the product in the table . 


'ELE ( elel ele2 - addr ) \ compute address of product 

ID - SWAP ID - GORD * + \ offset from base address 

GBLK BLOCK GOFFS + \ address of start of table 

+ ; \ address of element 
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Now we can define a word to print the ta¬ 
ble of a group. We first define a word that 
takes two elements and returns their 
product: 

: G* ( elel ele2 - elel*ele2 ) 

'ELE C@ ; 

Notice we have named this word g* sug¬ 
gesting it computes the product of two 
elements in a group. 

: .ELE ( ele — ) 

\ print element with a trailing space. 
EMIT SPACE ; 

The following word produces the limits for 
loops over all elements in the current 
group. 

: GLIMITS ( — u 1 ) 

\ loop limits to run through group 
GORD ID + ID ; 


: TABLE ( n — ) >GROUP CR 
2 SPACES GLIMITS 

DO . " _" LOOP CR 

GLIMITS DO I .ELE ." 

GLIMITS DO JIG* .ELE LOOP CR 
LOOP CR ; 

You can't usually tell 
much about a group 
by looking at its table. 
I have found that the 
ability to generate ta¬ 
bles in various for¬ 
mats is a useful way 
to transport informa¬ 
tion to other com¬ 
puter programs. 

Return to Top 


Example: 

8 table 


A |A B C D E F 
B |B C D E F A 
C |C D E F A B 
D |D E F A B C 
E |E F A B C D 
F |F A B C D E 
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Input and Output 


At this point we have a way of represent¬ 
ing and storing the main pieces of data for 
the program - the group tables. We can 
"install" a group as the current group. We 
have defined the main "group theory" 
word, g*, which determines multiplication 
in the current group. We have the output 
word table to display the table of a given 
group. In some sense, the word >group is 
an input word: It allows us to specify the 
number of the group currently of interest. 

The next step is to develop words that 
compute group information by chasing 
around the table. We need a way to input 
an element of the current group. We can¬ 


not, at this moment, type "B C G*" to get 
the product of B and C in the current 
group. We could, however, type 66 67 
G*. WHY? First note that G* takes two pa¬ 
rameters that are expected to be ASCII 
codes of upper case letters within the 
range of the current group. The parame¬ 
ters for G* must be on the stack in ad¬ 
vance. If we type B C G*, the system will 
start reading from the left and will try to 
find B and C as words in the dictionary. If 
it cannot find them, it will try to interpret 
them as numbers in the current base (10). 
So 66 67 G* will work - and the letters 
won't. 
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Here is what happens on the current Forth 
system (Win32Forth) running the original 
code: 

8 >group ok 
BCG* 

Error: C is undefined 

In this case, B is actually a word in the dic¬ 
tionary (having something to do with the 
editor) - but its action is not to put 66 in 
the stack! On the other hand, C is neither 
a word in the dictionary nor a recognized 
number. 4 


Keep in mind that G* is expecting ASCII 
codes for group letters. An input format 
that requires the user to remember and 
put in these codes is not friendly. There 
are several ways to handle the input of 
group letters. 


4 Forth systems print "ok" after completing execution of 
a line of input to let the user know that the commands 
have been processed and the system is waiting for input. 
In this article the "ok" prompts have been removed from 
the printout for most examples. 
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Input of Elements - method 1 (using a handler) 


Forth is perfectly happy to work with 
strings that are not names for words in the 
dictionary (and not numbers): We just 
precede them by a "handler," a word that 
reads a string from the following input and 
does something specific with it. 

The Forth word ASCII (which is replaced 
by CFIAR in the current Forth Standard) 
takes a following string and returns the 
ASCII code of its first character. In the 
following example we use the word . (dot) 
to print the number on top of the stack 
(and remove it). 


ASCI I A . 65 

ASCI I B . 66 

ASCI I b . 98 

ASCI I C . 67 

66 EM T B 
98 EM T b 

So, for example, we could say ASCII B 
ASCII C G* to multiply elements B and C. 
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In the 1990 version, I introduced a new 
handler, GETE: 

: GETE BL WORD 1+ C@ ; 

The Forth word WORD does not take pa¬ 
rameters from the stack - it reads the in¬ 
put stream. WORD, preceded by a delim¬ 
iter, takes the next string from the input 


stream up to the delimiter and returns the 
address as a counted string. 

So gete gets the next word (up to blank) 
and returns its first character. It is used 
just like ASCII but has a name that sug¬ 
gests that we are getting a group element. 
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Input of Elements - method 2 (making group elements dictionary words) 


If we wish to use the simple syntax B C 
G* we could just make A, B, C, etc. words 
in the dictionary whose action is to put 
their ASCII codes on the stack. The easi¬ 


est way to 
stants: 

65 CONSTANT 
67 CONSTANT 
69 CONSTANT 
71 CONSTANT 


do this is to make them con- 


A 

66 

CONSTANT 

B 


C 

68 

CONSTANT 

D 


E 

70 

CONSTANT 

F 


G 

72 

CONSTANT 

H 

etc 


There is a major disadvantage to this: The 
words I and J are the names used in Forth 
for loop indices. If they are redefined, they 
cannot be used with their original mean¬ 
ings. Any particular Forth implementation 
may already have some useful meaning 
attached to other letters. 5 


5 There is a vocabulary mechanism in Forth which makes it possi- 
ble to have different actions for the same name (it's is- it's like hav¬ 
ing a set of dictionaries). I actually tried this first. This mechanism 
is suitable if a new meaning is needed for a period of time (e.g. the 
assembler and editor in many systems are put in separate vocabular¬ 
ies). It is not convenient to frequently switch back and forth between 
vocabularies, so this proved to be an inconvenient solution. 


Return to Top 
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A better approach is to use slightly differ- This allows the syntax /B /C G* when 

ent names, /a, /b etc for the group ele- typing from the keyboard. 


ments: 


65 CONSTANT /A 

67 CONSTANT /C 

69 CONSTANT /E 

71 CONSTANT /G 

66 CONSTANT /B 

68 CONSTANT /D 

70 CONSTANT /F 

72 CONSTANT /H etc. 
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Some First Procedures 


: EORDER ( ele — ord ) ID GLIMITS 
DO OVER G* DUP ID = 

IF 2DROP I ID - 1+ LEAVE THEN 
LOOP ; 


This computes the order of an 
element by taking powers until 
the identity is found. 
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Output using EORDER for 
group 8 

Table of group 8 for reference 

definition of EORDER 

comments for definition 

/A EORDER . 1 
/B EORDER . 6 
/C EORDER . 3 
/D EORDER . 2 
/F EORDER . 6 

8 table 

: EORDER ID 

GLIMITS DO 

(x x A n ) 

OVER 

G* 

DUP ID = 

IF 

2DROP 

I ID - 1+ 

LEAVE 

THEN 

LOOP ; 

loop over ASCII codes 
<- loop invariant 
( x x A n x ) 

(x x A [n+1]) 

If pwr is ident 

compute order 

A |A B C D E F 
B|BCDEFA 

C ICDEFAB 

D |DEFABC 

E |EFABCD 

F |F A B C D E 


: POWERS GETE ID GLIMITS 

DO DUP .ELE OVER G* DUP ID = 
IF 2DROP LEAVE THEN LOOP ; 

This computes and prints the powers of an 
element in the given group. Note that the 
group letter comes after the powers com¬ 
mand: 


POWERS A A 
POWERS B A B C D E F 
POWERS C A C E 
POWERS DAD 
POWERS EAEC 
POWERS F A F E D C B 
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There are two groups of order 4 (numbers 4 and 5 in our list of groups). They are not iso¬ 
morphic - which can be immediately seen by computing the orders of their elements. 


4 >GROUP 

5 >GROUP 

/A EORDER . 1 

/A EORDER . 1 

/B EORDER . 2 

/B EORDER . 4 

/C EORDER . 2 

/C EORDER . 2 

/D EORDER . 2 

/D EORDER . 4 


At this point I intended the program just to 
find a quick way to determine whether two 
groups of order 1-16 are isomorphic. An 
obvious step toward solving the isomor¬ 
phism problem is to recognize that certain 
properties of a group (such as the order of 
the group and the order of its elements) 
must be preserved. It seemed useful, 
therefore, to write a command to loop 
through all elements in a group and print 
their orders. 


Here is a very simplified version of OR¬ 
DERS. In the next section we will develop 
fancier versions: 

: ORDERS ( grp# — ) >GROUP 

GLIMITS DO I EORDER . LOOP ; 

4 ORDERS 1222 

5 ORDERS 1424 

Return to Top 
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The Evolution of ORDERS 


In this section I discuss several attempts 
to improve the functionality of ORDERS. 
The process illustrated here, where several 
versions of a word are produced, each im¬ 
proving upon the previous in some re¬ 
spect, is common. We always have the 
older version available for comparison - 
the system is never "broken" at any point. 
We try to write so that individual words 
can be improved without affecting other 
parts of the system. 

The ORDERS command provides important 
information, but group letters may be re¬ 
arranged in an isomorphism, so the orders 
of elements may be in a different order for 


two isomorphic tables. It may not be easy 
to see at a glance if two groups have the 
same number of elements of each order. 
It would also be much easier to compare 
two groups if the outputs were in a match¬ 
ing form. One approach would be to sort 
the list of orders of elements. A simpler 
approach, since we just want the number 
of elements of each order, is to define a 
new version of ORDERS that does book¬ 
keeping. We need an array to save the 
number of times each order occurs. Here 
we look at a way to build special types of 
array structures by allocating memory and 
producing words to access it. 
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I wrote the original version of Groupsl6 
using F83, a 16-bit version of Forth (inte¬ 
gers occupied 16 bits or 2 bytes) so a se¬ 
quence of integers in memory had ad¬ 
dresses 2 units apart. To allocate a block 
of memory for storage of up to 20 integers 
I used 

CREATE OCNT 40 ALLOT 

Invoking the word ocnt puts the base ad¬ 
dress of this block of storage on the stack. 
We now need to locate the address of the 
k th integer in this block: 

: 'OCNT ( k - addr ) 2* OCNT + ; 

We need to be able to increment the num¬ 
ber at this location: 


: +OCNT ( k — ) 'OCNT 1 SWAP +! ; 

And we need to be able to initialize the ar¬ 
ray (setting all integers to zero): 

: OCLEAR OCNT 40 ERASE ; 

This code contains assumptions about the 
size of an integer and is one of the few 
places where modification of the code was 
needed when I transported it to 
Win32Forth for this article. Win32Forth is a 
32-bit version, and addresses for consecu¬ 
tive integers are 4 bytes apart. The "bit- 
tedness" assumptions in the earlier code 
are in the use of 40 ALLOT to make space 
for 20 integers and in the use of 2* in the 
calculation of the address for a particular 
table entry. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 52 


What is bad about this is that one cannot 
update the code automatically. (Some¬ 
times 2* is used for multiplication by 2 - it 
cannot be just changed to multiplication by 
4 everywhere.) 

Here is a way to make the code "imple¬ 
mentation independent": Rather than in¬ 
troduce "magic numbers" in the code, we 
introduce a word or words to contain the 
implementation-specific information used 
for allocation and address calculations. 
The word CELL is used to refer to the size 
of an integer or stack element. For F83 it 
has the value 2, and for Win32Forth it has 
the value 4. The word CELLS multiplies a 


number on the stack by CELL. 6 With this 
understanding, the following definitions 
would work on either version of Forth: 

CREATE OCNT 20 CELLS ALLOT 

: 'OCNT ( k - addr ) CELLS OCNT + ; 

: +OCNT ( k — ) 'OCNT 1 SWAP +! ; 

: OCLEAR OCNT 20 CELLS ERASE ; 

The use of cells rather than 2* (which 
does the same in F83) has the advantage 
of alerting the reader to the fact that 
memory allocation is the issue here, rather 
than arithmetic. It can be useful to have 
different names for the same action to 
clarify the purpose for which the action is 
being taken. 


6 These words are part of the ANS-Forth standard, so 
are already part of modern Forth systems and do not 
have to be defined by the user. 
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We now define an improved version of 
ORDERS: 

: ORDERS ( grp# — ) >GROUP 

OCLEAR 

GLIMITS DO I EORDER +OCNT LOOP 

GNUM 3 .R GORD 3 .R 2 SPACES 
20 1 DO 

GORD I MOD 0= 

IF I 'OCNT @ . THEN 

LOOP ; 

There are two sections of this code, one to 
count the orders and the second to print 
the results. Notice that we print the num¬ 
ber of elements of order k only when k is a 
divisor of the order of the group. 

8 orders 8 6 1122 

21 orders 21 12 138000 


A slight further refinement of the output 
section makes it possible to produce a 
comma- and quote-delimited file for use in 
a spreadsheet. For future reference, here 
are words that print a comma or quotation 
marks, respectively: 

: COMMA ASCII , EMIT ; 

: QUOTES ASCII " EMIT ; 

These can be used to produce "comma- 
delimited" output that is accepted by most 
spreadsheet programs. The spreadsheet 
file looks like this: 
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gnum 

gord 

# elements of 
each order 

1 

1 

ii ii 

2 

2 

" 1 1" 

3 

3 

" 1 2" 

4 

4 

" 1 3 0" 

5 

4 

" 1 1 2" 

6 

5 

" 1 4" 

7 

6 

" 1 3 2 0" 

8 

6 

"112 2" 

9 

7 

" 1 6" 

10 

8 

" 1 7 0 0" 

11 

8 

" 1 3 4 0" 

12 

8 

" 1 5 2 0" 

13 

8 

"116 0" 

14 

8 

"112 4" 


The ability to export these results to a 
spreadsheet program (where they can be 
sorted and examined) proved to be a ma¬ 
jor step in identifying which of the original 
set of groups might be isomorphic. The 
number of elements of each order is pre¬ 
served under isomorphism - groups that 
differ in this invariant are certainly not iso¬ 
morphic. We ask whether the following 
theorem is true: 






















JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 55 


Theorem? Two groups of order n are 
isomorphic if and only if they 
have the same number of ele¬ 
ments of each order. 

I think most of us would regard such a 
theorem as unlikely. (If it were true, it 
would have been proved by Lagrange or 
Cauchy, and we would find it in text¬ 
books.) If it is false, what is the first 
counter-example? 


Notice that the chart above shows that up 
to order 8, groups are distinguished by the 
numbers of elements of various orders. 
The theorem is also true if we require the 
groups to be abelian. (Two abelian groups 
are isomorphic if and only if they have the 
same number of elements of each order.) 
For what N is the theorem true for n < N? 

My last step in modifying ORDERS was to 
produce output in a more pleasant format. 


4 orders 

5 orders 

Group number 4 of Order 4 

Group number 5 of Order 4 

1 elements of order 1: A 

1 elements of order 1: A 

3 elements of order 2: BCD 

1 elements of order 2: C 

0 elements of order 4: 

2 elements of order 4: B D 
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Producing this output requires saving not 
only the count but also the order of each 
element. The original code, as mentioned, 
had built-in assumptions that integers 


were 16-bit. We show both the original 
and modified version of the code for this 
final version of ORDERS. 


Original version for F83 - 24Mar90jjw 

Revised for Win32Forth 

CREATE OTABLE 40 ALLOT 

CREATE OC NT 40 ALLOT 

CREATE OTABLE 20 CELLS ALLOT 

CREATE OC NT 20 CELLS ALLOT 

' ORD ( el e - - addr ) 

ID- 2* OTABLE -I- ; 

' ORD ( el e - - addr ) 

1 D - CELLS OTABLE -I- ; 

' OCNT ( i - - addr ) 

2* OC NT -I- ; 

' OCNT ( i - - addr ) 

CELLS OC NT -1- ; 

■+OCNT ( ord -- ) 

' OC NT 1 SVNAP -H ; 

TOC NT ( ord -- ) 

1 OC NT 1 SVNAP -H ; 

ORD! ( ord el e - - ) 

' ORD ! ; 

ORD! ( ord el e - - ) 

' ORD ! ; 

OCL EAR 

OTABLE 40 ERASE 

OCNT 40 ERASE ; 

OCL EAR 

OTABLE 20 CELLS ERASE 

OCNT 20 CELLS ERASE ; 
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Using these auxiliary words, the new defi¬ 
nition of ORDERS is: 

: ORDERS ( grp# — ) CR >GROUP 

." Group number " GNUM . 

." of Order " GORD . CR OCLEAR 
GLIMITS DO 

I EORDER DUP +OCNT I ORD! LOOP 
GORD 1+ 1 DO GORD I MOD 0= 

IF 3 SPACES I 'OCNT @ 2 .R 
elements of order " I 2 .R : " 

GLIMITS DO I 'ORD @ J = 

IF I .ELE THEN 
LOOP CR 
THEN LOOP ; 

This code is hard to read. Too many differ¬ 
ent functions are combined in one word. It 


takes a sharp eye to discern that part of 
this code is devoted to computing orders 
and part to printing results. Code that is 
designed to produce aligned output is 
mixed with code to perform computations. 
Code that is essentially routine "boiler¬ 
plate" conceals the main ideas. By the 
time of Groups32, I paid some attention to 
writing style, as I will discuss shortly. 

Here is the current version of ORDERS: 

: Orders ( grp# — ) CR >Group 

." Group number " Gnum . 

." of Order " Gord . CR 

Calc-Orders Prt-Orders ; 


Return to Too 
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Finding Isomorphisms 


All of the groups in Almquist's list have dif¬ 
ferent distributions for the orders of ele¬ 
ments until we get to his groups of order 


12. Here we find, for the first time, two 
groups in his list that have the same num¬ 
ber of elements of each order. 


Gnum 

Gord 

#elements of each < 

25 

12 

" 1 

1 

2 

2 

2 

4" 

24 

12 

" 1 

1 

2 

6 

2 

0" 

22 

12 

" 1 

3 

2 

0 

6 

0" 

20 

12 

" 1 

3 

8 

0 

0 

0" 

21 

12 

" 1 

3 

8 

0 

0 

0" 


Groups #20 and #21 are, therefore, the 
first groups in the original set that could 
be isomorphic. To show that they are iso¬ 
morphic, we need to display an isomor¬ 
phism between them. A first attempt at 
displaying an isomorphism is to create a 


mechanism for showing what happens to a 
table under a mapping of elements. We 
imagine a mapping X that sends the letters 
A, B, C, ... to the same letters in a per¬ 
muted order. 
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CREATE 'XG 20 ALLOT \ This will contain the images of A, 

: XG ( ele - xele ) \ given ele this returns the image xele 
ID - 'XG + C@ ; 

: XG' ( xele - ele ) \ this finds the inverse image of a letter 

GLIMITS DO I XG OVER = 

IF DROP I LEAVE THEN LOOP ; 

: XG* ( xl x2 - x3 ) \ Multiply in image group 

XG' SWAP XG' SWAP G* XG ; 

: >XG ( follow by string in letters ) 

\ this allows input of the mapping 
BL WORD COUNT 'XG SWAP CMOVE ; 


: RTABLE ( n — ) \ This prints the table of the image group 

>GROUP CR 

2 SPACES GLIMITS DO ." _" LOOP CR 

GLIMITS DO I .ELE 

GLIMITS DO J I XG* .ELE LOOP CR 
LOOP CR ; 
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>XG ACBD 

This applies map A^-A, 
B->C, C->B, D->D 

5 TABLE 

A |A B C D 

B | B C D A 

C | C D A B 

D |D A B C 

Original table for group 5 

_A_C_B_D 

A |A C B D 

C |C B D A 

B | B D A C 

D | D A C B 

Intermediate table ob¬ 
tained by interchanging 

B and C wherever they 
occur. 


5 RTABLE 

Final table obtained by 


rearranging rows and col¬ 

A |A B C D 

umns 

B |B A D C 


C | C D B A 


D |D C A B 



We have applied the isomorphism X and 
obtained the resulting table. Can we apply 
an isomorphism like this to Group #20 and 
transform its table to that of Group #21? 
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There are 12! = 479,001,600 candidates 
for X. However, if X is an isomorphism 
from group 20 to group 21, it must send A 
to A, it must send an element of order k to 
an element of order k, and it must send a 


product in the first group to the product of 
the images in the second group. (In fact, 
the map is determined by what it does to 
the generators - but we do not yet have 
tools to determine the generators.) 
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If we map B to B, then we must map C to 
C, since C is B*B in both groups. D must 
get mapped to an element of order 2. 
Let's try D-»E. B*D is F in Group 20, B*E 
is F in group 21 so we should map F-»F. 
Continue working in this fashion. 
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This can be seen to be the same as the ta¬ 
ble for group 21. We have, therefore, 
shown that groups 20 and 21 are isomorphic 
by displaying an isomorphism between 
them. 


The isomorphism command in Groups32 developed from this 
idea. It is much more sophisticated: It allows the user to 
experiment with images of elements (undoing trials if they 
do not seem to work), automatically computes and fills in all 
consequences of a partial mapping, checks for inconsisten¬ 
cies, etc. 


Return to Top 
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Comments on Improvements 


It is possible to criticize and improve upon 
almost everything done so far. While it is 
good to reevaluate and improve code, it 
has proved to be bad to do this prema¬ 
turely. The first goal should be to produce 
something that works and to experiment 
with options. A working system is the most 
important tool for developing improved 
versions. The modularity of Forth- 
produced systems often makes it possible 
to revise segments of a system while re¬ 
taining code that comes before and after. 

It is often not possible, or even desirable, 
to pre-conceive a system and design it on 


paper before writing any code. Forth does 
not require this (an advantage for this 
type of work). It usually becomes appar¬ 
ent only when a system is used what parts 
of it need speeding up, what parts could 
use a more congenial nomenclature, what 
functions are really needed, how they 
should work, and what form is best for the 
information they produce. This approach to 
system design can be called "middle out" 
coding: A rudimentary version of a system 
is produced, and it is refined and expanded 
as it is being used. 
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Changing the language to fit the program 

I preferred, for this article, to transport the 
earlier code to Win32Forth, which is used 
to run the current Groups32. The earlier 
Forth system used 16-bit integers and the 
current system uses 32-bit integers. In the 
earlier version, the group data were stored 
on disk using a virtual memory mecha¬ 
nism. In the modern version, the data are 
stored in memory as part of the dictionary. 
This change in the way data are stored did 
not actually require altering the code for 
the rest of the system - it was handled by 
an approach that may be unique to Forth: 
Instead of changing the code, change the 
language. We define a word BLOCK that 
reads data stored in a segment of mem¬ 


ory. It appears to the rest of the code to 
act just as if the older virtual memory 
scheme was being used: 

\ + pseudo block structure + 

CREATE PSEUDO 12 1024 * ALLOT 
: BLOCK ( n — addr ) 

1024 * PSEUDO + ; 

When the program is loaded, group data 
are read from the original block file into 
this segment of memory. The rest of the 
program does not know that the definition 
of BLOCK has been given a new meaning. 
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Choosing better names for words 

Now let me talk about some of the things I 
didn't like and that were changed. For one 
thing, I didn't like the names that I gave 
to some things in the original version. The 
word GLIMITS, for example, is used to put 
the loop limits on the stack for a loop over 
all elements of the current group. The 
etymology is obvious - but it is a name for 
how the word does something rather than 
what it does. In use I found myself always 
having to look up what name I had used - 
there is something not at all memorable 
about the name "GLIMITS". In Groups32, 
this has been replaced by the word 
"For-AII-Elements". For similar reasons, 


the word EORDER has been replaced by 
EleOrder. 7 


7 Some readers might like "ElementOrders" better. The 
choice of names is a matter of taste. Very terse names 
like "EO" tend to make code hard to read. Very long 
names can also be hard to read. I prefer "Orders", for ex¬ 
ample, to 

"Display_the_Number_of_Elements_of_Each_Order". 
The maximum length of a name in Win32Forth is 63 
characters - so this would be legal but unwise. 
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When I have used this code with others 
new to Forth, I found that it was helpful to 
distinguish words that come as part of the 
Forth language and those that have been 
added as part of the Groups package. 
Words that are part of Forth are typed in 
caps, words added for the system are in 
mixed case. Contrast: 


: ORDERS ( grp# — ) >GROUP 

GLIMITS DO I EORDER . LOOP ; 

: Orders ( grp# — ) >Group 

For-All-Elements 
DO I EleOrder . LOOP ; 

As I have used Forth, the importance of a 
good choice of names has become clear. 
With proper choices, one can often read a 
section of code without having to refresh 
one's memory by looking up a description 
of the word in a glossary or in the source 
code. A properly chosen name for a word 
can often suggest what it does. 
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Inconsistent methods of action 

I also do not like inconsistencies among 
actions of words. To enter elements as 
parameters, we have several possibilities. 
In the code above we have words that re¬ 
quire the ASCII code of group letters in 
advance - and we introduced special 
words such as /A, /B, /C to allow this to be 
done conveniently. The word EORDER 
uses this approach: We say /B EORDER 
to get the order of B. The word POWERS, 
however, is based on the handler GETE. 
We say POWERS B to get the powers of 
B. Both are words that take an element of 
the current group as a parameter. Both 
approaches are valid. 


I have found it confusing to have some 
words that take the parameter in advance, 
and some that take it afterwards. This 
practice requires that you remember this 
matter of placement for each word. It is 
better to have similar words act consis¬ 
tently, if possible. Words that act on ele¬ 
ments should all either take their argu¬ 
ment in advance or take it afterward. It is 
more consistent with the style of Forth to 
have parameters for the word occur in ad¬ 
vance - so, eventually, I eliminated the 
approach based on GETE. It is easy to 
change the definition of POWERS to take 
its argument from the stack: 
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Syntax: POWERS B 

Syntax: /B POWERS 

: POWERS GETE ID GLIMITS 

DO DUP .ELE OVER G* 

DUP ID = 

IF 2DROP LEAVE THEN 
LOOP ; 

: POWERS ID GLIMITS 

DO DUP .ELE OVER G* 

DUP ID = 

IF 2DROP LEAVE THEN 
LOOP ; 


When faced with this kind of choice, it is 
usually not wise to spend a lot of time in 
meditation before writing the code. The 
best approach is to implement both possi¬ 
bilities and use them for a while before de¬ 
ciding the winner. (In keeping with the 
evolutionary theme of this article, one 
adopts a "survival of the fittest" approach 
to modifications.) 
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Improving speed 

Notice that the group tables store letters 
A, B, C, ... for the table entries. Since 
finding the product of two elements in¬ 
volves calculating a location in memory, 
every use of G* involves first subtracting 
ID (the ASCII code for the identity) from 
the two letters. It should be possible to 
use integers 0, 1, 2, ... for the internal 
representations of the group elements, 
which would eliminate some extra calcula¬ 
tions for G*. It does require, however, 
that group elements be converted to and 
from ASCII upon input and output. 
Speeding up G* (used many times within a 
computation) is far more important than 
speeding up printout (which occurs only at 
the end). Input and output operations oc¬ 


cur relatively infrequently and tend to be 
sufficiently fast. 

I should also be mention at this point that 
one can define selected Forth directly in 
assembly language. The word G* is exe¬ 
cuted very frequently, since most words in 
the system use it either directly or indi¬ 
rectly. In the following chart, I compare 
execution times for three different defini¬ 
tions for G*: the original Groupl6 defini¬ 
tion (group elements stored as letters), a 
Groups32 version using a "high level" 
Forth definition, and the current Groups32 
version using an assembly language defini¬ 
tion. 
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The timings were obtained using the same due entirely to the way the data are repre¬ 
version of Forth (Win32Forth) running on sented and the way G* is coded , 

the same computer - so differences are 


Timing for 10,000,000 executions of G* 

Original Groupsl6 // table stores 
group letters A, B, C, ... 

98420 ms 

Current Groups32 ver 7.0 // table 
stores 0, 1, 2, ... 

35750 ms 

Same as above with G* written in 
assembly language 

7300 ms 
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Of course, the value of trying to improve 
the speed of a single word depends on 
how much gain can be achieved, how fre¬ 
quently the word is used, and whether or 
not the additional performance is needed. 
If a word is used mainly interactively (i.e. 
a command usually issued from the key¬ 
board) and it returns an answer in 1/10 
second, it is hardly worthwhile to try to get 
it to execute in 1/100 second. A word like 
G*, however, is executed quite frequent- 
ly. 8 


8 The most time-consuming command in Groups32 
searches all groups of orders 1-32 for those having given 
generators and relations. In an example using this com¬ 
mand which took about 1.5 minutes, the G* word was 
executed 176,777,081 times. 


One should check code to see if the same 
computation is made repeatedly and 
should be replaced by code that computes 
once and stores the result. A number can 
be retrieved from a table far more quickly 
than it can be found by looping through a 
list. An example of this is the word XG' 
(badly named) that finds the inverse of an 
element under a map. In order to compute 
the product XG*, the inverse must be 
found twice. In order to display the trans¬ 
lated table (RTABLE), the product XG* 
must be computed for each pair of ele¬ 
ments. It is possible to rewrite >XG (the 
word that stores the map) so that it com¬ 
putes and stores the inverses. I made im¬ 
provements like this in the eventual ISO¬ 
MORPHISM command in Groups32. 
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Forth Style 

Groups32 is a software system that is de¬ 
signed to be modified and expanded as it 
is used. We will gradually produce a collec¬ 
tion of "top level" commands to be used 
interactively. These are supported by a 
collection of foundational and auxiliary 
words. 

The fact that Forth programming results in 
a large collection of words enhances the 
ability of a user to extend and modify the 


system. Forth style typically emphasizes 
the production of words with fairly short 
definitions. Deeply nested loops and con¬ 
trol structures can be avoided by introduc¬ 
ing auxiliary words. By keeping definitions 
short and uncomplicated, it is easier to 
spot errors. Equally important for system 
development is the fact that words can be 
executed interactively - one can test indi¬ 
vidual words in piecewise fashion as they 
are being written. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 74 


This process of breaking definitions into smaller parts is called "factoring". The most 
complicated word we have written so far is ORDERS: 

: ORDERS ( grp# — ) CR >GROUP 

." Group number " GNUM . 

." of Order " GORD . CR OCLEAR 
GLIMITS DO I EORDER DUP +OCNT I ORD! LOOP 
GORD 1+ 1 DO GORD I MOD 0= 

IF 3 SPACES I 'OCNT @ 2 .R 

elements of order " I 2 .R : " 

GLIMITS DO I 'ORD @ J = 

IF I .ELE THEN 
LOOP CR 

THEN LOOP ; 
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Notice that some factoring has taken 
place: Clearing the tables (OCLEAR), 

computing the order of an element (EOR- 
DER), incrementing the count (+OCNT), 
storing the order (ORDI), and printing an 
element (.ELE) have all been delegated to 
separate words. When writing ORDERS, 
we need to be able to perform these tasks, 
but we can forget the details of how they 
were performed. 

It is certainly possible to factor ORDERS 
still further. Some Forth writers like to 
make the top-level word look like this: 

: ORDERS ( grp# — ) 

Install-Group 
Initialize-Tables 
Compute-Counts 
Display-Info ; 


There is another reason for factoring: to 
separate sections of code that may be in¬ 
dependently useful. Here is a hypothetical 
word that, for some reason, needs to de¬ 
termine if the current group is abelian. The 
highlighted code accomplishes this by test¬ 
ing whether each pair of elements com¬ 
mutes. 

: BigWord 

<some stuff> 

TRUE 

GLIMITS DO 
GLIMITS DO 

J I G* I J G* <> 

IF DROP FALSE THEN 
LOOP LOOP 

<more stuff> ; 
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Making this a separate word will clarify the 
definition of BigWord, and it will also pro¬ 
duce a word that is potentially useful else¬ 
where: 

: Abelian? ( — t/f ) 

\ is current group abelian? 

TRUE 

GLIMITS DO GLIMITS DO 
J I G* I J G* <> 

IF DROP FALSE THEN 

LOOP LOOP ; 


Then we replace the code segment in Big- 
Word by the word Abelian? 

: BigWord 

<some stuff> 

Abelian? 

<more stuff> ; 
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Eliminating Magic Numbers 

A feature of good style in any program¬ 
ming language is the elimination of "magic 
numbers," that is, numbers that have spe¬ 
cial significance to the system or program 
and that appear in the code as literal 
numbers. Examples of magic numbers are 
16 (which appears as the maximum order 
of a group in Groupsl6), or 49 (which is 
the number of tables in the original Alm- 
quist posting), or 20 (which is used to allot 
storage for arrays - probably because it 
was bigger than 16 and seemed like a nice 
number). 

It is preferable to define magic numbers as 
constants and use the names rather than 
the numbers. This has the effect of making 


code more readable. The reader (most of¬ 
ten the author) does not have to wonder 
why a certain number appears in the code. 
MaxOrder more clearly indicates that (in 
the current system) the maximum order of 
a group is being used. This is preferable 
to using 16, which could be the maximum 
order of a group - but, then again, some¬ 
thing else. Or that may be the maximum 
order now, but not when the system is ex¬ 
tended. The elimination of magic numbers 
is also helpful if the code is used in a dif¬ 
ferent setting. In the present case, the 
system was eventually changed to handle 
groups of order up to 32 and it was trans¬ 
ported from a Forth system with 16-bit 
stack elements to one with 32. 
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We have already discussed the use of the 
word CELLS. This concept (with a different 
name) was introduced to me by the Kitt 
Peak VAX-Forth. At one time Forth was 
used very heavily by astronomers, and Kitt 
Peak produced a version of Forth allowing 
code to be written to run, without change, 
on 16 and 32-bit systems. I avoided a lot 
of grief in converting Groupsl6 to 
Groups32 by the fact that magic numbers 
appeared mainly in a few low-level words. 


It is wisest to write code assuming that 
you will want to understand it yourself in 
10 years, and perhaps modify it and use it 
for a different purpose. Mathematicians 
tend to treat computer code like theorems: 
They prefer not to re-prove something 
they already have proved. 


Table of Contents 
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Back to the Isomorphism Problem 


We can distinguish groups of order up to 
15 by the orders of the elements. The only 
groups in Almquist's list with the same set 
of element orders are groups 20 and 21 
(of order 12), and we have found an expli¬ 
cit isomorphism between them. 

Groups of order 16 provide the first coun¬ 
terexample to the hypothetical theorem. 
Groups 43 and 45 have the same numbers 
of elements of each order, but group 43 is 
abelian, while group 45 is not. 

Further isomorphisms occur in Almquist's 
tables among groups of order 16. There 
are also groups of order 16 that are not 


isomorphic but that agree in element or¬ 
ders and in commutativity (abelian or not). 
Thus we still do not have enough proper¬ 
ties to classify groups of order 16 up to 
isomorphism. 

Here are some additional, easily com¬ 
puted, properties that may be useful - I 
wrote them as experiments to find useful 
invariants. They do not appear in either 
Groupsl6 or Groups32, since they did not 
seem to be of lasting interest. They were, 
in fact, used to distinguish the original set 
of group tables. 
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Number of Commuting Pairs 

Eventually (in Groups32) we compute the 
center of a group and can easily compute 
the size of the center. At this point we 
just run through all pairs of elements and 
count those pairs that commute: 

: #Commuting ( grp — n ) >GROUP 
0 

GLIMITS DO GLIMITS DO 
JIG* I J G* = 

IF 1+ THEN 
LOOP LOOP ; 

Number of Squares 

We can define #SOUARES to count how 
many elements in a group are squares: 


CREATE SQ 20 ALLOT 

: #SQUARES ( grp# — n ) >GROUP 
SQ 20 ERASE 
GLIMITS DO 

I I G* ID - SQ + 1 SWAP C! 
LOOP 

0 

20 0 DO SQ I + C@ + LOOP ; 

If one actually wishes to see which ele¬ 
ments are squares, the following word will 
display the squares: 

: SQUARES ( grp# — ) >GROUP 

SQ 20 ERASE 
GLIMITS DO 

I I G* ID - SQ + 1 SWAP C! 
LOOP 

20 0 DO SQ I + C@ 

IF I ID + .ELE THEN 
LOOP ; 
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The task that motivated this work was 
accomplished at this point. 


It turns out that, at this point, we have 
implemented enough properties to classify 
groups up to order 16. Coupled with the 
mechanism for exhibiting isomorphisms, 
we can detect which of Almquist's tables 
are for isomorphic groups. (There are ac¬ 


tually only 42 groups of orders 1-16 up to 
isomorphism.) 

The following tables were found to be iso¬ 
morphic: (20,21), (33,40), (36,38), 

(32,39,42), (34,41). In the process, I 
produced a method that will quickly iden¬ 
tify any group of order 1-16 up to isomor¬ 
phism. Table of Contents 
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Groups'! 6 - Transition 

From this point on, the numbering of groups corresponds to the interme¬ 
diate revisions of Groups'!6 and to Groups32. The duplicate groups have 
been eliminated, and the numbering for groups of orders 1-16 is 
changed. Eventually groups of orders 17-32 will be added. 


In the transition to Groups32, I made a 
number of changes. ID appears so often 
in the code above because we stored 
group elements as A, B, C, ... , and we of¬ 
ten need to use the corresponding num¬ 
bers 0, 1, 2, ... to access tables or arrays. 
We have already pointed out that the 
group multiplication word G* would be 
faster if this conversion from the letters to 
indices were eliminated. I changed the 
group tables to store 0, 1, 2, ... for the 
elements, which required revising many of 


the existing words that use tables and ar¬ 
rays. GLIMITS (which put loop limits on 
the stack) was renamed to For-AII-Ele- 
ments, and the limits are now indices 
rather than ASCII codes. The words /A, 
/B, /C, ... that put ASCII codes on the 
stack were changed to put 0, 1, 2, ... on 
the stack. I changed the display word .ELE 


from 

.ELE 

( ele — ) 

EMIT 

SPACE ; 

to 

.ELE 

( ele — ) 

ID + 

EMIT SPACE ; 
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Reminder: Code in red was in the 
early 1990 version, Groupsl6; 
Code in green is in the current 
Groups32; Blue represents code 
that was not in the 1990 version 
and did not survive to Groups32. 


These changes were not as disruptive as it 
might seem at first. Many of the words in 
the system are not aware of how the data 
are represented and so required no 
change. 

A consideration in implementing a system 
is precisely this: Can the programming be 
done in such a way that details of data 
representation are hidden in lower level 
words? Do most of the words in the sys¬ 
tem need to know whether group elements 


are represented by indices or by ASCII 
codes? Do most of the words in the sys¬ 
tem need to know where or how the data 
are stored? 

These questions have repercussions when 
changes are made in the system. They 
also have psychological repercussions: 
How much is the mathematician allowed to 
think in terms of mathematical objects 
(groups, elements, etc.), and how much 
must he or she think in terms of computer 
representation. The purpose of some of 
the changes in this transition is to make 
the language and concepts more in terms 
of mathematical objects and procedures 
applied to them. Even the small change in 

name from GLIMITS to For-All-Elements is 

significant. 
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I had an important change in attitude in 
the transition: Several procedures were 
originally written to produce output on the 
screen - they were regarded as top-level. 
I found while extending the system that it 
would be useful to separate the generation 
of information from printing it. An example 
of this is a factoring of the Orders com¬ 
mand (discussed previously) into a proce¬ 
dure to calculate the orders of elements 
and a procedure to print them: 

: Orders ( grp# — ) CR >Group 

." Group number " Gnum . 

." of Order " Gord . CR 

Calc-Orders Prt-Orders ; 


This allows other commands to use infor¬ 
mation about orders of elements. This, in 
turn, required the development of data 
structures that could be used to transmit 
information from one procedure to an¬ 
other. In the case of subgroups and sets, it 
involved representing data in a form that 
could be transmitted on the stack. In the 
case of element orders, it meant providing 
other procedures access to the locations in 
memory where data were stored. 

The transition also was marked by the ad¬ 
dition of other features - most notably, a 
way to deal with subgroups and a permu¬ 
tation group package - plus refinement of 
existing features. 
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First Subgroups Package 


My subgroups code represents a first effort 
to introduce some of the further structure 
of group theory into the system. It is 
transitional code and does not represent 
the way that subgroups are handled in 
Groups32. It was, however, very useful, 
providing a way to test later code as it was 
being developed. It displays the subgroup 
generated by a collection of elements. In 
the transition period, I had not yet intro¬ 


duced the idea of treating a set as an ob¬ 
ject that could be put on the stack. 

In this package, a particular subgroup (the 
current subgroup) is treated as a counted 
table or array. This approach was eventu¬ 
ally replaced by a more sophisticated one. 
While the code in this section does not ap¬ 
pear in Groups32, it does mark an impor¬ 
tant stage in the transition. 


Overview 

Top Level words and Examples 

Details Return to Table of Contents 
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Overview 

The main purpose of this code is to find a 
way to exhibit the subgroup generated by 
a collection of elements. The subgroup 
being generated is represented by a 
counted array or table of bytes. This is a 
list of the elements found (so far) to be in 
the subgroup. The first byte is the count, 
and the remaining bytes are subgroup 
elements. 

The most important top level words in the 
package are SGfirst ( ele -- ) and SGnext 
( ele — ), which are used to enter the first 
and subsequent elements. The main 
words used in implementation are 
>SGtable ( ele -- ) and SGsweep. 

>SGtable adds a new element to the table 


(subgroup) if it is not already there. 
SGWEEP repeatedly computes products of 
pairs of elements currently in the sub¬ 
group, and uses >SGtable to add the prod¬ 
ucts to the subgroup. This continues until 
no new elements are added. 

Return to Top 
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Top Level words and Examples 


Top Level Words 

: SGfirst ( ele — ) 

SG-init >SGtable 

SGsweep SG. ; 

Use this for first ele¬ 
ment to initialize the 
table 

: SGnext ( ele — ) 

>SGtable SGsweep SG. ; 

Use this for remaining 
elements 


When I first introduced these words, I was using the conventions of the 1990 version. In 
the 1990 system, they had a slightly different form, using GETE to get the input element 
from the input stream - so one would use "SGfirst B" to start a subgroup with element B. 
The examples below were produced using the current version of 
Groups32. It is necessary to add the constants. 

0 CONSTANT /A 1 CONSTANT /B 

2 CONSTANT /C 3 CONSTANT /D 
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4 

CONSTANT 

/E 

5 

CONSTANT 

/F 

6 

CONSTANT 

/G 

7 

CONSTANT 

/H 

8 

CONSTANT 

/I 

9 

CONSTANT 

/J 

10 

CONSTANT 

/K 

11 

CONSTANT 

/L 

12 

CONSTANT 

/M 

13 

CONSTANT 

/N 

14 

CONSTANT 

/o 

15 

CONSTANT 

/P 

16 

CONSTANT 

/Q 

17 

CONSTANT 

/R 

18 

CONSTANT 

/s 

19 

CONSTANT 

/T 


The input is in the form "/B SGfirst". This output also assumes the new numbering of 
groups. 

32 >GROUP 

/B SGfirst B C D A 
/E SGnext BCDAEFGH 

/I SGnext BCDAEFGHI JKLMNOP 

The subgroup generated by B, E, and I is the whole group. Return to Top 
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Details 


We are storing the elements of a subgroup 
in SGtable. The first byte of SGtable is the 
number of elements already in the table. 
The others are the stored subgroup ele¬ 
ments (the data) in whatever order they 


have been generated. Let us say we have 
already found that elements B, F, and E 
are in the subgroup and have been found 
in that order. The SGtable will look like 
this: 


3 

B 

F 

E 


* 1 
addrl ac 

dr end 


The word SGtable puts the address of the addr and return addr_end (the next avail- 

table on the stack (marked addr). The able position at the end of the table), 

standard Forth word COUNT takes the ad¬ 
dress of a counted structure and returns 
the pair address_of_data and count. Thus 
COUNT would return the pair addrl 3 in 
this case. The sequence COUNT + will take 
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We would like to add a new element to the 
table if it is not already there. We use a 
clever (but well known) trick to do this. 
The strategy is to put the new element in 
the position at the end of the table and 
then scan to see if it is found. It will be 
found, of course. 

I am torn, in this section, between clarity 
and authenticity. I have decided to pre¬ 
sent code in the form in which it appears 
in some of the transitional versions of the 
system. I have, however, presented the 
code with component sections separated 
by blank lines. For clarity, >SGtable should 
have been factored. The structure is 


: >SGtable ( ele — ) 

\ add ele to table if new 

DUP 

Put_at_end 

Find_position ( pos ) 

Found_at_end? ( t/f) 

IF Update-count THEN ; 

There are obvious inefficiencies in the 
code: SGsweep, for example, recomputes 
products that have already been com¬ 
puted. It really should take only products 
with any newly added elements. Remem¬ 
ber, however, this code is intended to dis¬ 
play a subgroup as elements are added 
from the keyboard - so its speed is really 
limited by the input and output. 
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Auxiliary Words 

CREATE SGtable 40 ALLOT 


: SG. SGtable COUNT 0 

?DO DUP C@ .ELE 1+ LOOP 

DROP CR ; 

Print the subgroup 

: SG@ ( k - ele [k] ) 

SGtable 1+ + C@ ; 

Return the k-th element, k=0, 1, 2, ... . 

: SG-init SGtable 40 ERASE 

} 


: >SGtable ( ele — ) 

DUP SGtable COUNT + C! 


Put the element at the end of the table. 

SGtable COUNT + 1+ 

SGtable 1+ 

DO ( ele ) 

I C@ OVER = 

IF DROP I LEAVE THEN 

LOOP 


Loop over data addresses. 

Check if stored element is ele. 

If so, put address where found on stack. 

SGtable COUNT + = 

IF SGtable C@ 1+ 

SGtable C! THEN ; 


Check if found at end of table (where we 
stored it) - if so, increment count of table. 










JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 92 


: SGord ( — sgord) 

SGtable C@ ; 


: SG* ( kl k2 - prod ) 

SG@ SWAP SG@ SWAP G* ; 


: SGsweep 

SGord 0= 

ABORT" Table is Empty " 

BEGIN SGord 

SGord 0 DO SGord 0 DO 

J I SG* >SGtable 

LOOP LOOP 

SGord = UNTIL ; 

Take products of existing elements in 
SGtable and add to table. Repeat this 
until no new elements are added. 

Top Level Words 

: SGfirst ( ele — ) 

SG-init >SGtable 

SGsweep SG. } 

Use this for first element to initialize 
the table. 

: SGnext ( ele — ) 

>SGtable SGsweep SG. ; 

Use this for remaining elements. 


Return to Top 
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Sets 


In previous sections, the top level words 
have printed information to the screen. 
The information they produced was not in¬ 
tended to be used by further words. In 

Overview 

Top Level Words and Examples 

Details 

Sets and Subgroups 


this section we introduce a representation 
for sets (subsets of group elements) as ob¬ 
jects that can be transmitted to other 
words on the stack. 


Return to Table of Contents 
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Overview 

In October, 1991 I wrote a new version of 
Groupsl6, using a newer and more exten¬ 
sive 16-bit version of Forth (FPC). This re¬ 
write included some of the changes above. 
The most notable new idea, however, was 
to introduce tools for representing and 
manipulating sets. Some of the commands 
could be written to produce a set as a re¬ 
sult, rather than print a list of elements. 
Sets could be transmitted on the stack 
from one procedure to another. 


Sets are represented by "bit maps": stack 
elements (the same as integers) in which a 
bit is set when an element is in the subset. 
For example, the set {B, D} contains ele¬ 
ments 1 and 3 (A is element 0) and so is 
represented by 10, which is 1010 in bi¬ 
nary. This allows the representation of 
subsets of {A, B, .... ,P}. Procedures are 
provided for input and output, as well as 
standard set manipulations, test of set 
membership, etc. 


Return to Top 
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Top Level Words and Examples 


#Set 

( set — tele) 

Returns number of elements in 

set 

2 A n 

( n — 2 A n ) 

Raise 2 to the input power. 

(See Singleton). 

Contains 

( setl set2 — f ) 

Returns true if setl contains 
set2 . 

Empty-Set 

( — 0 ) 

Returns the empty set 

Member! 

( ele set — set' ) 

Make ele a member of set 

Member? 

( ele set — f ) 

Returns true if ele is a member 

of set 

Set- 

( subl sub2 — subl-sub2 ) 

Returns difference of sets 

Set & 

( subl sub2 — intersection ) 

Returns intersection of sets 

Set* 

( setl set2 — product ) 

Returns the set of products of 
elements in setl set2 in given 
group. 
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Set. 

( set — ) 

Prints set (given as integer 
bitmap) using letters for ele¬ 
ments 

Set + 

( subl sub2 — union ) 

Returns union of sets 

Set-G 

( — gset ) 

Returns underlying set of cur¬ 
rent group 

Singleton 

( n — {n} ) 

Returns set with ele as only 
member 

SubGrp. 

( subg — ;;; print subg ) 

Prints subgroup (bitmap as 
set) using letters. This is 
the same as Set. 
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Examples 

{ ABC} { BC} Contains . -1 (true) 

{ ABC} { BD} Contains . 0 (false) 

/B { ABC} Member? . -1 

/D { ABC} Member? . 0 

/D { ABC} Member! 

Set. { A B C D } 

/B Singleton Set. { B } 

{ ABC} { BCD} Set- Set. { A } 

{ ABC} { BCD} Set+ Set. { A B C D } 

{ ABC} { BCD} Set& Set. { B C } 

{ ABCD} #Set . 4 


Return to Top 


8 TABLE 

_A_B_C_D_E_F_ 

A |A B C D E F 
B |B C A F D E 
C |C A B E F D 
D |D E F A B C 
E |E F D C A B 
F |F D E B C A 
Group number 8 of Order 6 

1 elements of order 1: A 

3 elements of order 2: D E F 

2 elements of order 3: B C 

0 elements of order 6: 

{ B} Generates Set. {ABC} 

{ BD} Generates Set. {ABCDEF } 
Center Set . { A } 


/B Centralizer Set. {ABC} 

{ B} Generates Normalizer Set. {ABCDEF } 
{ D} Generates Normalizer Set. {AD} 
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Details 

: 2 A N ( n - 2 A n ) \ integer with bit n set 

1 SWAP 0 ?D0 2 * LOOP ; 

: Singleton ( n - set ) 2 A N ; 

: Contains ( setl set2 - t/f ) \ TRUE if setl contains set2 
SWAP OVER AND = ; 

: Member? ( ele subg - t/f ) 

SWAP Singleton Contains ; 

: Member! ( ele subg — subg 1 ) 

SWAP Singleton OR ; 

Notice the use of "bitwise" operations AND and OR. 5 2 OR gives 7 while 5 2 AND gives 0. 

5 (10 1) 

2 (010) 
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Here is a word for printing a set: 

: SubGrp. ( subg — ) \ print subg 
." { " 1 ( bit mask ) 

16 0 DO 2DUP AND \ check if bit is set 
IF I .ELE THEN 
2* \ left shift mask 

LOOP 2DROP ." } " ; 


Notice that the idea of representing sets 
(and subgroups) as arrays has been aban¬ 
doned in favor of using a representation as 
integers. The word SubGrp. here is de¬ 


fined to print a subgroup represented as 
an integer - it supercedes the early ver¬ 
sion of SG. in which a set is represented as 
an array - but it has the same function. 


: Set. SubGrp. ; 
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Here is a word for inputting a set. The 
definition uses some features of Forth 
compilation to make it "state smart" (i.e., 
it acts differently if the system is compil¬ 
ing). Notice that, when interpreting, this 
word produces an integer, which is put on 
the stack. When the system is compiling, it 


stores the integer in the definition. This 
definition is not included in Groups32, 
since the command completion interface 
used in Groups32 prompts for input of 
sets. The details of this definition are be¬ 
yond the scope of this article. 


: { ( — set ;;; interpret string up to } as a set ) 

0 \ start with empty set 

ASCII } WORD \ get string up to } 

COUNT 2DUP UPPER \ convert to upper case 

0 DO DUP C@ ID - DUP 0 15 BETWEEN \ check if valid 
IF ROT Member! SWAP \ set addr 

ELSE DROP THEN 1+ \ next address 

LOOP DROP 

STATE @ IF [COMPILE] LITERAL THEN ; IMMEDIATE 
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Some basic set operations: 

: Set- ( setl set2 — setl-set2 ) 

\ set difference 

OVER AND XOR ; 

: Set+ ( setl set2 — union ) 

\ set union 
OR ; 

: Set& ( setl set2 — intersection ) 
\ set intersection 
AND ; 

: #Set ( set — tele ) 

0 SWAP 

BEGIN ( cnt n ) DUP 
WHILE 0 2 ( cnt lsb msb 2 ) 

UM/MOD ( cnt r n/2 ) 

>R + R> \ update cnt 
REPEAT DROP ; 


In #Set we shift bits to the right by divid¬ 
ing by 2. The Forth word 

UM/MOD ( ud un - ur uq ) 

takes a double precision unsigned integer 
and divides by a single precision integer. It 
produces an unsigned remainder and un¬ 
signed quotient. In each step in the loop 
we are dividing by 2, the remainder (0 or 
1) is added to the count, and the quotient 
is used in the next step. 

A double precision integer is represented 
by two stack elements; the top of the 
stack is the most significant part. 
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An important operation on subsets of a 
group is to form S*T = { s*t | seS, teT }. 

Here is the high-level Forth version: 

VARIABLE Stemp 

: Stemp! ( ele — ) Stemp @ Member! Stemp ! ; 

: Set* ( setl set2 — setl*set2 ) 0 Stemp ! 

For-All-Elements DO OVER I SWAP Member? 

IF ( first ele is member of setl ) 

For-All-Elements DO I OVER Member? 
IF ( second ele is member of set2 ) 
JIG* Stemp! 

THEN LOOP 

THEN LOOP 2DR0P Stemp @ ; 


Return to Top 
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Sets and Subgroups 


We are now in a position to write a word 
"Generates" that finds the subgroup gen¬ 
erated by a given set. Notice that the idea 
here is the same as used previously in 
SGsweep: We repeatedly take all products 
of elements in a subgroup under construc¬ 
tion until we get something closed under 
products. In this case we take only prod¬ 
ucts using new elements. The subgroup is 
represented as a set, so adding new ele¬ 


ments uses the set operations we have 
just defined. We use the following notation 
in the comments for Generates: 

S is the originally given set of generators. 

X is the total set of products formed so 
far. 

Y is the set of products yet to be multi¬ 
plied by S. 
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: Generates ( set — subg ) >R 


1 1 ( 

BEGIN R@ Set* ( 

OVER Set- ( 

SWAP OVER Set+ ( 

SWAP ?DUP ( 

0= UNTIL 

R> DROP ; 


X = Y = {id} ) 

X Y*S ) 

X Y' Y'=Y*S-X new prods ) 

Y' X' X'=X+new ) 

leave if no new products ) 


Since Set-G, the underlying set of G, occurs in subgroup operations, we have >Group cal¬ 
culate and save it when a new group is installed. 
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: LCoset ( ele subg — lcos ) SWAP Singleton SWAP Set* ; 

: RCoset ( ele subg — rcos ) SWAP Singleton Set* ; 

: Centralizer ( ele — centrz ) 1 

For-All-Elements 

DO OVER I OVER G* SWAP I G* = 

IF I SWAP Member! THEN 
LOOP SWAP DROP ; 

: Center ( — center ) 1 

For-All-Elements 
DO I Centralizer Set-G = 

IF I SWAP Member! THEN LOOP ; 

: Normalizer ( subg — normal ) 1 

For-All-Elements 

DO OVER I OVER LCoset I ROT RCoset = 

IF I SWAP Member! THEN LOOP SWAP DROP ; 

: Normal? ( subg - t/f ) 

Normalizer Set-G = ; 


Return to Top 


Return to Table of Contents 
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Listing Subgroups 


We wish to make a list of all distinct sub¬ 
groups and also maintain a list of the cor¬ 
responding generators. Thus the ordered 
list will be indices into a table of subgroups 
and generators, subgtable is a record 
structure that holds subgroups and gen- 

Ordered Lists Package 

We make use of a package for handing or¬ 
dered lists that was written in 1989. 
Rather than exhibit the code for this pack¬ 
age, we describe the words used here. 


erators for them. The comparison function 
sgcomp was chosen so that subgroups get 
listed in increasing size. After generating 
all subgroups, one can access them as an 
array subgroup with indices 

0..#SUBGROUPS-l . 


Note: The implementation of ordlist used here uses 
single bytes for indices and therefore can contain at most 
256 indices. 
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ORDLIST 

<number> ORDLIST <name><comparison> 

This is a defining word. 
<number> is the maximum 
number of elements in the 
list. <comparison> is the 
name of a comparison func¬ 
tion with the stack diagram 
(i j —1 | 0 | 1 ) 

>LST 

(list - ) 

Install list as the current list 

ISIN? 

( n list - i t/f) 

Check if the element of index 
n is in the list. If it is, i is the 
position, and the flag is 

TRUE. If it is not on the list, 
then i is the place where it 
should be inserted. 

(INS) 

(ni~) 

Insert n at position i 

LST@ 

(i-n) 

n is the entry at position i in 
the current list. 

OL-EMPTY 

(list - ) 

Empty the given list. 


An ORDLIST is paired with an 
array-like structure that 
holds the actual data. The 
ordered list is a list of indices 
of data. When we create an 
ordered list, we supply a 
comparison function. 

Suppose i and j are two indi¬ 
ces and L[i], L[j] the corre¬ 
sponding data. The compari¬ 
son function has stack dia¬ 
gram 

(ij~ -1 I 0 | 1). 

-1 if L[i] < L[j], 

0 if L[i] = L[j], and 
1 if L[i] > L[j]. 

One of the uses of ordered 
lists in this system is to keep 
a list without duplicates. 
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CREATE SUBTABLE 100 2* CELLS ALLOT 9 
: 'SUBG ( index — addr_of_subg ) 

2 CELLS * SUBTABLE + ; 

: 'GEN ( index — addr_of_gen ) 'SUBG CELL + ; 

: SGCOMP ( idxl idx2 — —1|0|1 ) 

'SUBG @ SWAP 'SUBG @ SWAP 
OVER #SET OVER #SET - ?DUP 0= 

IF 2DUP U< IF 2DROP 1 ELSE U> THEN 

ELSE >R 2DROP R> 0> IF -1 ELSE 1 THEN 
THEN ; 

100 ORDLIST SUBGLIST SGCOMP 


In the following code, we establish 
an array SUBGTABLE that holds 
subgroups and generators. The 
comparison function SGCOMP com¬ 
pares subgroups. SI < S2 if SI has 
fewer elements than S2 or, if they 
have the same number of elements, 
SI comes before S2 as unsigned in¬ 
tegers. SUBGLIST is the name for 
the ordered list. When this name is 
used, an address is put on the 
stack. >LST uses this address to 
make SUBGLIST the current ordered 
list and all the ordered list words re¬ 
fer to it. 


VARIABLE SGPTR 

: >SGTABLE ( genset — ) \ put at end 

DUP 1 SGPTR +! SGPTR @ 'GEN ! 
GENERATES SGPTR @ 'SUBG ! ; 


9 In Groups32 this array was made somewhat larger to handle 
groups up to order 32. 
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: NEWSUB ( ele subgidx — 
2DUP 'SUBG @ MEMBER? 


NEWSUB takes an element and the index of a 
subgroup in SUBTABLE. It will produce a new 
subgroup if the element is not in the existing one. 


IF 2DR0P 

ELSE 'GEN @ MEMBER! 

>SGTABLE 

SGPTR @ SUBGLIST ISIN? 


IF DROP -1 SGPTR +! 

ELSE SGPTR @ SWAP 
(INS) THEN 


if ele is in the subgroup, then drop both 

Else add the element to generating set. This returns a 
new generating set with the element added. 

>SGTABLE takes a generating set and puts it and the 
subgroup it generates at the end of SUBTABLE. 

We now check the ordered list SUBGLIST to see if this 
subgroup has previously appeared. 

If it is already there, we DROP the position (produced 
by ISIN?) and decrease the SGPTR - in effect, 
removing the subgroup from the end of the table. 

If it is new, we insert the index in the proper place in 
the ordered list SUBGLIST. 


THEN 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik- 110 


GENERATE-SUBGROUPS repeatedly tries to 
add generators to existing subgroups. No¬ 
tice that the inner loop is over all ele¬ 
ments, while the outer loop is over only 


the indices of SUBTABLE corresponding to 
subgroups added in the previous pass. 
(Any subgroups prior to that already have 
been processed.) 


: GENERATE-SUBGROUPS ( grp# — ) 

>GROUP SUBGLIST OL-EMPTY 
SUBGLIST >LST 0 SGPTR ! 

1 0 'SUBG ! 0 0 'GEN ! 0 0 (INS) 

-1 

BEGIN ( old_top ) \ previous top of SUBTABLE 

1+ SGPTR @ DUP >R \ save current SGPTR value 
1+ SWAP \ use only recently added subgroups 

DO \ loop over subgroups 
FOR-ALL-ELEMENTS 
DO \ loop over elements 
I J LST@ NEWSUB LOOP 

LOOP 

R> SGPTR @ OVER = \ compare old SGPTR with current 

UNTIL ( exit when no more subgroups have been added ) 

DROP ; 
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Once the subgroups have been generated, 
we display them in increasing order, first 
ordered by size then by lexicographic order 
of elements. 

: SHOW-SUBGROUPS 
SUBGLIST >LST 

CR . 11 * = Normal subgroup" 

CR Generators" 25 TAB ." Subgroup" 

SGPTR @1+0 

DO CR I 3 .R 2 SPACES I LST@ DUP 
'GEN @ SET. 25 TAB 'SUBG @ 

DUP NORMAL? IF 42 ELSE BL THEN EMIT SG. 
LOOP CR ; 

After GENERATE-SUBGROUPS, the sub¬ 
groups are in an array. <k> 'SUBG gives 
the address of the subgroup in position k, 
0 <= k < N (where N is the number of 
subgroups). The index of the k-th sub¬ 
group in the ordered list is <k> LST@ 


: #SUBGROUPS SGPTR @ 1+ ; 

: SUBGROUP ( indx — subg or empty ) 

SUBGLIST >LST DUP 0 #SUBGROUPS WITHIN 
IF LST@ 'SUBG @ 


ELSE DROP 0 
: SUBGROUPS ( grp# — ) 

Example: 

8 SUBGROUPS 

* = Normal subgroup 
Generators 
0 { } 

1 { D } 

2 { E } 

3 { F } 

4 { B } 

5 { B D } 


THEN ; 

GENERATE-SUBGROUPS 
SHOW-SUBGROUPS ; 


Subgroup 
*{ A } 

{AD} 

{ A E } 

{ A F } 

*{ A B C } 

*{ A B C D E F } 


Return to Table of Contents 




JOMA 


Journal of Online Mathematics and its Applications 


Wavrik- 112 


Some Simple Programs 

Notice that we are developing a group the¬ 
ory language that can be used when we do 
further programming. Here are some ex¬ 
amples of new words created using this 
language: 

: Abelian? ( — t/f ) 

\ is current group abelian? 

CENTER SET-G = ; 

: IsAnOrder? ( n - t/f ) 

\ is n the order of some element? 

\ Assume CALC-ORDERS has been run 
FALSE ( n flag ) 

MaxOrd 0 DO OVER I 'ORD @ = 

IF DROP TRUE THEN 
LOOP SWAP DROP ; 


: Cyclic? ( grp# — t/f ) 
>GROUP CALC-ORDERS 
GORD IsAnOrder? ; 
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When put in a loop over all groups, Cyclic? 
allows us to identify the cyclic groups. 
Here are some variants on the IsAnOrder? 
command: 


: OfOrder ( n - ele t | f ) 

\ find element of order n 
FALSE ( n flag ) 

For-All-Elements 
DO OVER I EleOrder = 

IF 2DROP I TRUE LEAVE THEN 

LOOP 

DUP 0= IF SWAP DROP THEN ; 

0 CONSTANT Empty-Set 
: OrderN ( n — set ) 

\ set of all elements of order n 
Empty-Set ( n set ) 
For-All-Elements DO 
OVER I EleOrder = 

IF I SWAP Member! THEN 
LOOP SWAP DROP ; 
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The following is an incorrect attempt to 
identify the dihedral groups. It looks for 
groups of order n that contain an element 
of order 2 and an element of order n/2 and 
that are non-abelian: 


A dihedral group certainly satisfies these 
conditions, but some groups satisfy these 
conditions without being dihedral. 


: Oops-Dihedral? ( grp# — t/f ) 
>GROUP CALC-ORDERS 
Abelian? NOT 
2 IsAnOrder? AND 
GORD 2 / IsAnOrder? AND ; 
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Here, for example, are the two non-abelian groups of order 8. One of them actually is the 
dihedral group, and the other is the group of quaternionic units. It should be easy to spot 
group 13 as the dihedral group - it has lots of reflections (elements of order 2) and only two 
rotations (elements of order 4). However both groups are non-abelian and have elements of 
order 2 and n/2. 


Group number 13 of Order 8 

1 elements of order 1: A 

5 elements of order 2: C E F G H 

2 elements of order 4: B D 
0 elements of order 8: 


Group number 14 of Order 8 
1 elements of order 1: A 
1 elements of order 2: C 
6 elements of order 4: B D E F G H 
0 elements of order 8: 
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An easy way to locate the dihedral groups 
in the current Groups32 system is to use 
the search command. A dihedral group has 
two generators x and y, x is of order 2, the 
order of y can be left unspecified, and the 
relation is xy = y'x where ' denotes in¬ 
verse. It is an interesting exercise at this 
point to try to write a word DIHEDRAL? 
with the stack diagram ( grp# -- t/f ) that 
determines if the given group is dihedral. 
Here is a solution: 

: Dihedral? ( grp# — t/f ) 

>GROUP FALSE 
For-All-Elements DO 
For-All-Elements DO 
I EleOrder 2 = 

J EleOrder Gord 2/ = AND 
J I G* J G* I = AND 

IF DROP TRUE THEN 
LOOP LOOP ; 


I will use this as an example to explore 
some issues in the performance of code. It 
is obvious that this code is not efficient in 
several ways: (1) It runs through all pairs 
of elements of the group - it should exit 
the double loop if it finds a pair that satis¬ 
fies the conditions. (2) If the first ele¬ 
ment, I, is not of order 2, it nevertheless 
chooses all possible second elements, J. 
(3) There are three conditions to be 
tested, and they are tested in sequence - 
but all three are evaluated (and combined 
using AND), even if the first is false. The 
code, however, is direct, simple, and easily 
written. It has actually been optimized for 
clarity rather than for execution speed. 
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Be aware of the modern perversion that 
everything must run as fast as possible. 
Your time is valuable! This is a word that 
will identify, interactively, only one group 
for each even order - and it will do so in 
less than one second. It is simply not 
worth the programming time required to 
make this particular word run faster. It 
may, nevertheless, be instructive to use 
this as an example for examining ways to 
increase speed. 


We introduce a way to test the conjunction 
of conditions that stops testing when one 
is found to be false. Here is a way to syn¬ 
thesize this behavior from existing control 
structures. 

We want to create words ?IF and ?THEN so 
that ?IF cccc ?THEN behaves as follows. 
We assume a flag on the stack. If this flag 
is TRUE, it is removed from the stack and 
cccc is executed (leaving behind a new 
flag). If this flag is FALSE, then cccc is not 
executed, and the FALSE flag is left on the 
stack. This is logically equivalent to the 
sequence IF cccc ELSE FALSE THEN. 
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: ?IF POSTPONE 10 IF ; IMMEDIATE 

: ?THEN POSTPONE ELSE 


This version of Dihedral? will now execute 
a condition only if the previous conditions 
have evaluated to TRUE. With this change 
it takes only 234 milliseconds to find the 
dihedral groups of order up to 32 rather 
than 411 milliseconds. 


POSTPONE FALSE 

POSTPONE THEN ; IMMEDIATE 


: Dihedral? ( grp# — t/f ) 


>GROUP FALSE 


For-All-Elements DO 


For-All-Elements DO 
I EleOrder 2 = 

?IF J EleOrder Gord 2/ 
?IF J I G* J G* I = 


?THEN 

?THEN 


IF DROP TRUE THEN 
LOOP LOOP ; 


10 The phrase POSTPONE xxx is used in the definition of 
an immediate word WWW to prevent the action of xxx 
from being compiled into WWW. Instead it causes WWW 
to compile the action of xxx into any word in which WWW 
is used. 
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We have previously used the word LEAVE 
to exit prematurely from a DO .. LOOP 
construct. LEAVE, however, exits only from 
the current loop. Forth provides a word 
EXIT that exits from the current procedure 
and a word UNLOOP that discards the cur¬ 
rent loop parameters. Thus 


: Dihedral? ( grp# — t/f ) 

>GROUP FALSE 
For-All-Elements DO 
For-All-Elements DO 
I EleOrder 2 = 

?IF J EleOrder Gord 2/ = ?THEN 
?IF J I G* J G* I = ?THEN 

IF DROP TRUE 

UNLOOP UNLOOP EXIT 11 
THEN 

LOOP LOOP ; 


11 It should be noted that EXIT exits not just from the 
control structure, it exits from the entire procedure. You 
cannot use EXIT to leave a loop and continue with later 
code in the same procedure. 
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Since we have gone this far, here is a ver¬ 
sion of Dihedral? that takes care of all 
three of the objections: 

: Dihedral? ( grp# — t/f ) 

>GROUP FALSE 
For-All-Elements DO 
I EleOrder 2 = 

IF For-All-Elements DO 
I EleOrder Gord 2/ = 

? IF I J G* I G* 

J = ?THEN 
IF DROP TRUE 

UNLOOP UNLOOP EXIT 
THEN 
LOOP 

THEN 
LOOP ; 


I suspect that this can be written and un¬ 
derstood only by someone who first writes 
one of the simpler versions. It finds all the 
dihedral groups up to order 32 in 213 mil¬ 
liseconds. 
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Now, if one frequently needs to determine 
whether a group is dihedral, there is a 
rather interesting further possibility: Gen¬ 
erate results by a method that might be 
slow - and then use the stored results 
rather than re-calculating. 

CREATE DIHEDS MaxTables ALLOT 
\ number of groups 

: 'DIHED ( grp# — addr ) 

DIHEDS + ; 

: Old-Dihedral? Dihedral? ; 

: Dihedral? ( grp# — t/f ) 

'DIHED C@ ; 


: Calc-Dihedral 

DIHEDS MaxTables ERASE 
For-All-Groups DO 

I Old-Dihedral? I 'DIHED C! 
Loop ; 

Calc-Dihedral 

Now the time for executing Dihedral? is 
reduced to .03 milliseconds, not counting 
the time it takes to create the table. The 
moral of the story is that one should first 
optimize for clarity and go further only if it 
is necessary. The Forth language does 
support a variety of ways to make things 
faster and also perhaps clearer, ranging 
from the ability to introduce new control 
constructs to the ability to store previously 
computed information. 


Return to Table of Contents 
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Groups32 


At this point, if you are using Windows, 
you should run the Groups32 program in a 
separate window. This will make it possible 
to see how the program works and to try 
examples. Examples can be copied and 
pasted from the article or a text editor to 
Groups32. If you want to do more ex¬ 
tended programming, Win32Forth uses 
plain text files that can be created by an 
external editor and loaded into Groups32 
by using INCLUDE <filename>. When 


Groups32 starts, it is in the commands 
completion interface (to be described 
later). To do programming, one must QUIT 
this interface. 

A useful Forth word is "see". SEE 
<name>, where <name> is a word in the 
dictionary, will decompile the dictionary 
entry for the word, allowing you to see the 
definition of the word. 


Reminder: Code in red is found in the early 1990 version, Groupsl6; Code in 
green is found in the current Groups32; blue is used to represent 
code that was not in the early version and did not survive to 
Groups32. 
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Permutations 


This package provides the basic represen- the group generated by a given set of 
tation and manipulation of permutations. permutations. 

It also provides a procedure for computing 

Overview 

Top Level Words and Examples 

Details 

Permutation Groups 
Bugs, Error Traps, Etc. 


Return to Table of Contents 
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Overview 

A permutation is represented as an array 
of bytes. If P is a permutation of l...n, 
byte 0 is used to store n and bytes l...n 
store P(l)...P(n). A permutation is repre¬ 
sented on the stack by the base address of 
this block of memory. 

Access words are provided to manipulate 
permutations. 

If P sends i -> j, then i P Sends returns 
j. To make P send i ->■ k, do k i P Send!. 

(The order is chosen to be consistent with Forth "k addr !" to 
store k at addr.) 


In the earliest version of this package, I 
used a defining word to declare a permuta¬ 
tion and set aside storage for it. I provided 
words for input, output, and multiplication 
of permutations. Eventually I added a 
mechanism to deal with the storage of in¬ 
termediate results to allow permutations to 
be used in chain computations. I also 
added a procedure to find the group gen¬ 
erated by a set of permutations. 

Return to Top 
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Top Level Words and Examples 


We start with the earliest attempt. Notice 
that not all of these words are in the cur¬ 
rent version of Groups32. We make a de¬ 
fining word MakePerm that allocates a block 
of storage for a named permutation. We 
will access the storage as an array of 
bytes. If the permutation is understood to 
permute 1, ... , n, we store n in byte 0. In 
general, byte k will contain the image of k 
under the permutation. The word pinit 
( n perm ) initializes a permutation to be 
the identity. The word p. ( perm ) prints a 
permutation. To store a permutation, I in¬ 
troduced the word p! The stack diagram 
is ( Xi .. x n n perm ), where x k is the image 
of k under the permutation. An initial ver¬ 


sion of a multiplication word, here called 
op* ( perml perm2 - result ), multiplies 
two permutations and puts the product in 
a named permutation called Result. 
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Here is how this looks in use: 

MakePerm PI \ Declare PI and P2 
MakePerm P2 \ as permutations 

12434 PI P! \ Store data 
2 3 4 1 4 P2 P! 

PIP. { 1243} \ print permutations 

P2P. { 2341} 


Reminder: 

Code in red is found in the 
Groupsl6; Code in green is found 
in the current Groups32; blue 
represents code that was not in 
the early version and did not sur¬ 
vive to Groups32. 


PI P2 oP* P. { 2314} 

\ use default left->right 
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We eventually add input and output in cy¬ 
cle notation, but the most serious problem 
at the moment is that the definition of 
multiplication puts all products in the same 
place: a permutation named Result. The 
problem can be seen in this example 

1 2 4 3 4 PI P! 

2 3 4 1 4 P2 P! 

PI P2 oP* DUP oP* P. { 3134} 

The correct answer is 

PI P2 P* DUP P* P. { 3124} 


What is happening is that the first product 
PI P2 oP* is stored as the permutation Re¬ 
sult. The sequence DUP oP* makes Re¬ 
sult both the input and output. Since oP* 
alters the permutation Result, it is alter¬ 
ing its input argument in this case. The 
conclusion is that for chain calculations 
there must be more than one place for re¬ 
sults. 
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It is important to understand that this is 
not a problem peculiar to permutations - it 
arises even for integer arithmetic. The 
stack provides the kind of flexible storage 
needed. If we perform a*b + c*d (using 
algebraic notation) we would do 
a b * c d * + (using reverse Polish nota¬ 


tion). The result of the first multiplication 
must be temporarily stored until we per¬ 
form the second multiplication. The stack 
diagrams for this sequence of operations 
show a*b remaining on the stack until the 
final addition. 
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The Temporary Storage Problem 

There are several possibilities for dealing 
with the problem of intermediate storage. 
The first, and worst, is to try to put the 
data for a permutation on the Forth pa¬ 
rameter stack. This approach would make 
a permutation occupy several stack cells. 
It has the effect of burying useful data and 
inhibiting the use of the standard stack 
operations. 

A second, more promising, approach is to 
create a special permutation stack. Since 
the same problem of intermediate storage 
arises with other types of algebraic ob¬ 
jects, and since a system might have sev¬ 
eral types of objects, this would require 
the user to introduce a new set of stack 


operations for each stack, and to be aware 
of the contents of several stacks while 
programming. This approach is psychologi¬ 
cally burdensome. 

In Groups32 and other algebra systems I 
have produced, I have represented objects 
on the regular Forth stack by an address. 
The address occupies a single cell - so it is 
manipulated by the usual stack operations. 
The user can think of a permutation being 
on the stack. SWAP can be used to ex¬ 
change two objects: two permutations, a 
permutation and an integer, etc. Thus, 
representing an object by its address al¬ 
lows a consistent way to integrate new 
types of data into the system. However, 
we still need to have a pool of locations for 
intermediate results. 
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An improved multiplication word, P*, will 
remove two (addresses of ) permutations 
from the stack and multiply them. It will 
then put the result in the next available 
temporary storage location and return the 
address of this location to the stack. 

In some languages, a very large pool of 
temporary locations is available - and 
computation proceeds, using these loca¬ 
tions, until the pool is exhausted. The sys¬ 
tem then attempts to identify storage loca¬ 
tions that are no longer in active use. Un¬ 
referenced locations are called "garbage" 
and the process of reclaiming them is 
called "garbage collection". There can be a 


noticeable pause in execution while the 
system collects garbage. 

We use a mechanism with storage for just 
16 temporary objects. It does a very 
small, quick garbage collection when these 
are exhausted. Each type of object has its 
own pool of 16 storage locations. A 
mechanism is provided for marking an ad¬ 
dress as available. When a computation is 
performed on a certain type of object, the 
result is placed in the next available stor¬ 
age location, and the address is returned 
to the stack. If no location is available, the 
system tries to find which locations are no 
longer in active use. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik- 131 


Cycle Notation 

The output word c. uses the same naive 
algorithm that one uses for computing 
products by hand. We keep track of which 
numbers have already been output. Start 
with all numbers unmarked, and mark any 
element that is sent to itself. Then start 
with the first unmarked element. Print an 
opening parenthesis, print (and mark) the 
elements in its orbit, and print a closing 
parenthesis. 

The input word (which I have called 
Perm: ) reads from the rest of the com¬ 
mand line. It is based on a word that reads 
a cycle (numbers between two parenthe¬ 


ses) and it multiplies the cycles it obtains. 
The cycles are not required to be disjoint. 
Perm: does check the syntax - so it re¬ 
ports errors such as improperly matched 
parentheses. 

Perm: (1 2) (1 3) 

C. (123) 

Perm: ((12) (13 

Syntax Error — Try Again 

? (1 2) (1 3) 

C. (123) 

Return to Top 
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Details 

33 CONSTANT PMaxSize 

\ largest n for permutations +1 for 
count 

: MakePERM CREATE PmaxSize ALLOT ; 

: 'Element ( i perm — addr ) + ; 

: Send! ( k i perm — ) 

'Element C! ; 

: Sends ( i perm — perm[i] ) 
'Element C@ ; 

: Degree ( perm — n ) 

0 SWAP Sends; 

: Degree! ( n perm — ) 

0 SWAP Send! ; 


: PInit ( n perm — ) 

2DUP 0 SWAP Send! SWAP 1+ 1 
?D0 I DUP 2 PICK Send! LOOP 
DROP ; 

: P. ( perm — ) 

\ print showing images of l..n 
. " { " DUP Degree 1+ 1 

DO I OVER Sends . LOOP 
DROP ." }" CR ; 

: P! ( xl .. xn n perm — ) 

2DUP Degree! 

1 ROT DO I OVER 

>R SWAP Send! R> 

-1 +LOOP DROP ; 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 133 


DEFER Direction 
MakePERM Result 

: oP* ( perml perm2 — result ) 

Direction 
DUP Degree 

DUP Result Degree! 1+ 1 
DO ( perm2 perml ) 

I 2 PICK Sends ( perm2 perml perm2[I] ) 
OVER Sends ( perm2 perml perml [perm2 [ I ] ] ) 

I Result Send! 

LOOP 2DR0P Result ; 

: Left->Right [ 1 ] NOOP IS Direction ; 

: Right->Left [ 1 ] SWAP IS Direction ; 
Left->Right \ default 

Recall that the multiplication operation de¬ 
fined here puts the product in a fixed per¬ 
mutation called Result and returns the ad¬ 


dress of that location to the stack. This 
code allows for the fact that some books 
use left to right multiplication while others 
use right to left. Forth allows for "vectored 
execution" using the word DEFER. The line 
"DEFER Direction" makes a dictionary en¬ 
try for Direction, which can then be used 
in subsequent code. An action for Direc¬ 
tion must be filled in before it is used. The 
code is written so that multiplication would 
be left to right without Direction (so we fill 
in the action NOOP when this is the desi¬ 
red direction). If if multiplication from 
right to left is desired, the action for Direc¬ 
tion is SWAP. The words Left->Right and 
Right->Left allow the direction to be swit¬ 
ched without recompiling. 
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Temporary Storage 

The user must understand that 
these temporary locations are a 
scratchpad area for operations on 
the objects - any result of lasting 
use must be moved to a perma¬ 
nent location. The garbage collec¬ 
tor assumes that an address is be¬ 
ing actively used if and only if it is 
on the stack. The internal work¬ 
ings of the temporary storage 
module are transparent to the 
user. Here is how the mechanism 
is introduced into the permuta¬ 
tions package: 


Temp-Storage Perm-Temps 

Establish a new Temp- 
Storage structure called 
Perm-Temps. 



CREATE Perm-Storage 

PMaxSize CELLS 16 * 

ALLOT 

Make a block of memory to 
store 16 permutations. 



: Perm-Setup 

Perm-Storage Perm-Temps 

16 0 DO DUP I Address ! 

PMaxSize CELLS + 
LOOP DROP ; 

Initialize Perm-Temps to con¬ 
tain the 16 addresses in the 
storage pool. 



Perm-Setup 

Execute the setup. 



: PTemp ( -- perm ) 

Perm-Temps Temp ; 

This is the main word that the 
user employs - it returns the 
address for the next available 
slot in the storage pool. 
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With this storage mechanism, it is very 
easy to fix P* without drastically changing 
the code. We need only to supply result 
with the address of the next available loca¬ 
tion in the storage pool. Notice that result 
may be used in several places in the code 
for an operation, and we do not want to 
use a different storage location at each 
reference. We therefore make result a 
VALUE (rather than a permutation). It will 
return the address of the storage area for 
a permutation. We define a word GetRes to 
store the address of the next temporary 
storage slot in result. 


0 VALUE Result 

: GetRes PTemp TO Result ; 

: P* ( perml perm2 — result ) 

Direction 

GetRes DUP Degree 

DUP Result Degree! 1+ 1 

DO ( perm2 perml ) 

I 2 PICK Sends ( perm2 perml perm2[I]) 
OVER Sends 

( perm2 perml perml[perm2[I]] ) 

I Result Send! 

LOOP 2DR0P Result ; 


Return to Top 




JOMA 


Journal of Online Mathematics and its Applications 


Wavrik- 136 


Permutation Groups 

The main remaining interesting feature of 
the permutations package is the ability to 
find the group generated by a given set of 
permutations. The algorithm for doing this 
is to maintain an array of the generating 
permutations (Gtabie) and an array of 
permutations generated (ptabie). We use 
the ordered list module, discussed previ¬ 
ously, to make the latter array an ordered 
list. Start with just the identity permuta¬ 
tion in the ptabie. Repeatedly multiply all 
elements in the ptabie by the generators - 


adding any permutation to the Ptabie that 
is not already there. Continue doing this 
until a pass (multiplying by generators) 
does not produce anything new. 

The following example shows how this fea¬ 
ture is used from the command line. I 
should mention that I eventually equip 
Groups32 with a user interface that makes 
it easier and quicker to use. 
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\ Klein's 4-group 

4 Size! \ This initializes the tables 
Gen: (1 2)\ add the generators 
Gen: (3 4) 

Make-Group\ apply the algorithm 
Show-Table\ a multiplication table 
\ can be obtained 


A B C D 
B A D C 
C D A B 
D C B A 

Show-Elements 

Group is of order 4 
A () B (3 4 ) C (1 2 ) 
D (1 2 ) (3 4 ) 
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At the moment, the permutations package 
seems disjoint from the rest of Groups32. 

We cannot immediately apply any of the 
procedures for orders, subgroups, etc. We 
can, however, generate a table. We now 

: InstallO \ Install as 

GSize @ Gord! 

GSize @ 0 
DO GSize @ 0 

DO J I PG* JIG*! LOOP 
LOOP ; 

If the order of the generated permutation 
group is between 1 and 32, we automati¬ 
cally have the table stored as table 0. (The 
slot for group #0 is used for temporary 
storage.) 


provide a mechanism for installing this ta¬ 
ble into the Groups32 system - replacing 
any of the tables 1-5 (which are for very 
small groups and will hardly be missed). 


table 0 


Groups32 provides Install ( n — ), which 
will move the table from slot 0 to replace 
table n, where n is 1, 2, 3, 4, or 5. Once a 
table is installed, all the commands can be 
applied to it. 
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Assume the Klein 

InstallO 

1 Install 
1 Table 

_a_b_c_d_ 

A |AB C D 
B |B AD C 
C |C D A B 
D |D C B A 


4 group is generated as above. 


\ move the current 
\ permgroup to slot 0 
\ install as group 1 
\ show the table 


1 Subgroups 

* = Normal 
Generators 
0 { } 

1 { B } 

2 { C } 

3 { D } 

4 { B C } 


subgroups 

Subgroup 

* { A } 

*{ A B } 

*{ A C } 

* { A D } 

*{ A B C D } 


\ show the 

subgroup 


Since Show-Elements provides a concordance of 
letters A, B, C, ... with permutations, it is easy to 
interpret any results in terms of permutations. 


Return to Top 
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Bugs, Error Traps, Etc. 

This is an appropriate place to discuss 
these issues, since I eventually equipped 
the permutations package with error han¬ 
dling at the point when I first made it 
available for use by others. I added the er¬ 
ror trapping even before I equipped the 
system with the user interface that I dis¬ 
cuss in the next section. The permutation 
package was not part of the original 
Groupsl6 - it was developed separately 
and eventually merged into the group the¬ 
ory package. 

When a system of this type is used by the 
person who wrote and implemented it, it is 
usually possible to minimize "bugs" that 
arise from coding errors (or from dis¬ 


agreements between the programmer and 
computer about how the world works). Let 
us call these "bugs of the first kind". Forth 
style emphasizes writing short, uncompli¬ 
cated definitions. Words can be tested in¬ 
teractively as soon as they are written. It 
can be quickly determined if a word does 
not function as intended, and the mistake 
is usually quickly pinpointed. 

When a piece of software is used by oth¬ 
ers, however, another type of bug can ap¬ 
pear: improper behavior resulting from in¬ 
correct use of the system. We will call 
these "bugs of the second kind". (Some 
call these "features".) 
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A perfect example of a "bug of the second Groupsl6. Recall the definition of g*: 
kind" occurs in the earliest version of 


\ compute address of product 
\ offset from base address 
\ address of start of table 
\ address of element 

: G* ( elel ele2 - elel*ele2 ) 

'ELE C@ ; 


: 'ELE ( elel ele2 - addr ) 
ID - SWAP ID - GORD * + 
GBLK BLOCK GOFFS + 

+ ; 


(ID is a constant, 65, the ASCII code of 
upper case A.) Group elements are repre¬ 
sented in this version by upper case let¬ 
ters. The code assumes that the elements 


are entered as ASCII codes of upper case 
letters - within the range of the particular 
group. 
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4 Table 

: Try GETE GETE G* ." = " EMIT ; 

A |A B C D 

Try B C = D 

B | B A D C 

Try b c = A 

C |C D A B 

Try E F = C 

D |D C B A 



When we put in B and C as uppercase, we 
get the correct answer. This word, how¬ 
ever, is case sensitive. When we put in b 
and c as lowercase, we do not get the cor¬ 
rect answer. We also get an answer when 
we put in E and F, even though they do 
not represent letters in the range of this 
group. 


There is no "bug of the first kind" in the 
definition of g*. When g* is supplied with 
valid input it does provide the correct out¬ 
put. However, the code for g* contains 
some assumptions about how it is to be 
used - and it contains no error traps to 
ensure it is used correctly. The errors seen 
here arise because g* was used incorrectly. 
This is a "bug of the second kind". 
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When the author is the user, he or she 
knows the intended form of input. It can 
be difficult for the author to anticipate 
some of the things other users might do - 
or how the software will respond. In the 
user interface (next section), I have tried 
to restrict "free form" input so that 
Groups32 can be used correctly by those 
unfamiliar with the system. The input for 
most commands is prompted, and it is 
easy to test if the input is of the proper 
type. For example, if the user must enter 
the number of a group, it is easy to check 
that the input is a number between 1 and 
144. If the user is to put in group ele¬ 
ments, it is easy to check whether they 
are valid. (We automatically convert to 
upper case - making the input case insen¬ 


sitive - but we do check for range.) These 
are cases where it is rather simple and 
quick to insure valid input. Many bugs of 
the second kind can be eliminated by 
prompting for input in a restricted form 
and rejecting invalid input. 

In some cases, however, we want users to 
provide input as a line of text - an exam¬ 
ple is the word Perm: in the permutations 
package. This word is supposed to be 
supplied with a string that represents a 
collection of cycles in the numbers l...n, 
where n has been specified in advance. 

Left->Right 
3 TO PSize 
Perm: (1 2)(1 3) 

C. (123) 
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Here is an early version of code for Perm: (the names of support words have been 
changed to correspond to the current ones). 


: 'Element ( i perm — addr ) + ; 

: Send! ( k i perm — ) 'Element C! ; 

: Sends ( i perm — perm[i] ) 

'Element C@ ; 

: Degree ( perm — n ) 0 SWAP Sends ; 

: Degree! ( n perm — ) 0 SWAP Send! 

: PInit ( n perm — ) 2DUP 0 SWAP Send! 
SWAP 1+ 1 

DO I DUP 2 PICK Send! LOOP DROP ; 

VARIABLE SvDepth 0 VALUE Result 
: GetRes PTemp TO Result ; 

: P* ( perml perm2 — result ) 

GetRes DUP Degree 

DUP Result Degree! 1+ 1 

DO I 2 PICK Sends OVER Sends 


I Result Send! 

LOOP 2DR0P Result ; 

5 VALUE PSize 

: << GetRes PSize Result PInit 
DEPTH SvDepth ! ; 

: >> DUP >R DEPTH SvDepth @ - 1 
DO OVER Result Send! LOOP 
R> Result Send! Result ; 

: Perm: 

\ use single parentheses in following string 
PTemp PSize OVER PInit 
BEGIN << 

ASCII ( PARSE 2DR0P 
ASCII ) PARSE EVALUATE 
DEPTH SvDepth @ - WHILE 
>> P* REPEAT ; 
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The size of a permutation (the n in Sn) is 
assumed to be in the VALUE PSize. The 
words << and >> are used for input of a 
single cycle. Here is what happens when 
we type << 1 2 3 >>: 

The word « saves the current depth of the stack and ini¬ 
tializes a new temporary permutation (called Result). The 
numbers 1 2 3 are just put on the stack. The word » com¬ 
pares the stored stack depth to the new stack depth, thus de¬ 
termining how many numbers are on the stack. A copy of 
the top number on the stack is saved on the return stack and 
a loop is entered that makes each number the image (under 
the permutation) of the number below it: 2 —> 3, 1 —> 2. 

The loop is completed by making 3 —> 1. This is how a cy¬ 
cle is handled. 


more cycles). It uses the << .. >> pair to 
interpret each cycle, and it multiplies the 
cycles it receives. Perm: ends its work 
when the input line is exhausted. 


The word Perm: reads the input that fol¬ 
lows it (presumably containing one or 
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This code functions correctly as described 
- as long as the user understands the as¬ 
sumptions built into Perm: 

(1) There must be a space between 
Perm: and the rest of the input line. 

(2) Perm: must be followed on the input 
line by a collection of cycles. 

(3) There must not be anything else on 
the input line. 

(4) Each cycle must begin with an open¬ 
ing parenthesis "(" and terminated 
with a closing parenthesis ")". 


(5) Within each cycle there must be 
numbers in the range l...n. 

(6) PSize must be set to n before Perm: 
is used. 

(7) The numbers within a cycle must be 
separated by spaces. 

(8) The numbers within a cycle must be 
distinct. 

(9) No characters, other than digits and 
spaces, can occur between the paren¬ 
theses. 
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Now, what happens if a user supplies 
Perm: with illegal input? Suppose the user 
does not close parentheses, does not 
match parentheses, or puts illegal charac¬ 
ters within parentheses? The system 
could react in several ways, ranging from 
aborting the command to producing some¬ 
thing that is not a permutation, and that is 
handled improperly by P*, C., and other 
words in the package. 

In the first example above, with g*, invalid 
input produced incorrect answers because 
numbers are fetched from the wrong part 
of memory. No damage is caused, just in¬ 
correct results. Some words, however, 
store things in memory, and it is possible 
that incorrect input could actually alter the 
code in memory. When Perm: is executed, 


the word Send! is ultimately called. This 
word has the stack diagram ( k i perm — ). 
It will store k in the memory position of 
index i in the given permutation. If i is not 
in the proper range, a number will be 
stored in memory other than within the al¬ 
located space for the permutation. If k is 
not within range, the result will not be a 
valid permutation, and subsequent opera¬ 
tions may lead to numbers being stored in 
improper places. Since Forth implementa¬ 
tions usually allocate data storage within 
the dictionary, it could turn out that num¬ 
bers can be erroneously stored in places 
being used by code - and the code will, 
therefore, stop operating correctly. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 148 


This should not be taken as an indication 
of fragility. In most cases, words that 
store things receive their arguments from 
other words that produce only valid input 
for them. 

The final versions of Perm: and its auxiliary 
words enforce the requirements for a cor¬ 
rectly formed input string. The input string 
is checked for balanced parentheses and 
for only digits within the parentheses. A 
range-checked version of Send! is used to 


make sure the numbers in a cycle are 
within range and that they are distinct. 
Within the command completion interface 
(see the next section), Perm: is used in a 
prompted input sequence, so the user 
must specify n first. (The word that stores 
n also checks that it is within range.) All of 
this is necessary just to try to eliminate 
bugs of the second kind. 
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Range Checking 

Any operation that stores data in memory 
could, if given erroneous input, store data 
in unintended places. To guard against 
this, one can add range checking to words 
that store data. Here is an example from 
code we have already examined. 

In the word Orders there is a word +Ocnt 
that increments the number in a certain 
slot in an array. Here is the original ver¬ 
sion: 

CREATE OCNT 33 CELLS ALLOT 
: 'OCNT ( k - addr ) CELLS OCNT + ; 

: +OCNT ( k — ) 'OCNT 1 SWAP +! ; 

Notice that k must be in the range 0-32 for 
this to work correctly. Now here is a ver¬ 
sion that checks that k is in range: 


CREATE OCNT 33 CELLS ALLOT 
: 'OCNT ( k - addr ) CELLS OCNT + ; 

: +OCNT ( k — ) 

DUP 0 32 BETWEEN NOT 

IF CR . ." Out of Range " ABORT 

ELSE 'OCNT 1 SWAP +! THEN ; 

There is, of course, a speed penalty for 
this extra checking. The original +OCNT 
executes 100,000 times in 26 milliseconds, 
while the range-checked version takes 50 
milliseconds. The ranged checked version 
is twice as slow. 
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Should range-checking be included? 

Range checking is obviously not necessary 
if a word is supplied its parameters only by 
other words that are guaranteed to supply 
valid parameters. +OCNT is a perfect ex¬ 
ample. This is a word that should NOT in¬ 
clude range checking. 

To see this, we need to look at how 
+OCNT is used. +OCNT occurs in the code 
for Orders, which calculates the orders of 


elements for one of the existing groups (of 
order 1-32). Thus, in use, the word 
+OCNT cannot receive an out-of-range in¬ 
dex. If we were to use the range-checked 
version, the number k would be checked 
each time to see if it is between 0 and 32, 
but the test will always be true because 
+OCNT will never receive a k that does not 
satisfy this condition. 
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One of the virtues of producing software 
for your own use is that it is not so neces¬ 
sary to deal with bugs of the second kind 
(and with error trapping). It can be very 
time-consuming to prepare a piece of 
software for use by others. Those who 
have lived with toddlers will appreciate 
that it is very much like "child proofing" a 
house. How much you have to do depends 
both on the complexity of your house and 
the sophistication of your toddler. Rather 
than put all your possessions in locked 
cabinets, however, you must give thought 
to what precautions are really most effec¬ 
tive. 


The critical words to consider for range 
checking (or other error trapping) are 
words such as Perm: that get input from 
the user and store things in memory. 
These are the words that have the poten¬ 
tial to alter the code itself. In the final ver¬ 
sion of the permutations package, I not 
only use a range-checked version of Send! 
but I also include other tests to make sure 
that the input string has the correct form. 
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At the beginning of this section, we used 
the example of G*. It would not be effi¬ 
cient to put range checking into a word 
like this, which is used very frequently. It 
would be better to examine the code that 
passes information to such words, making 
sure that G* will not receive invalid input. 
Once the internal code of the system is 
correct, the only source of invalid data will 
be input routines. This is usually the best 

Return to Top 


place to put error traps. Since data entry 
is slow anyhow, error traps here will not 
significantly slow the system. Once we are 
sure that the input routines will not send 
incorrect data to the rest of the system, it 
should not be necessary to have error 
traps on lower level words. 
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Search 


Groups32 has a mechanism to search the 
groups of orders 1-32 to find those with 
given generators satisfying given relations, 
or to find groups that contain subgroups 
with the given finite presentation. The al¬ 
gorithm used is quite literal: All group 
elements are substituted for the genera¬ 
tors, and then the orders and relations are 
checked. A "shortcut" evaluation is used: 
Once a condition is found to be FALSE, no 
subsequent conditions are checked. 

This is the slowest procedure in Groups32. 
Direct assembly language coding of G* 
(and 8 other words) does result in a dra¬ 
matic speedup of SEARCH. 


Here is how SEARCH can be used to iden¬ 
tify the dihedral groups, which are groups 
with two generators, x and y, with x of or¬ 
der 2, the order of y is not specified. The 
defining relation is xy = y'x, where ' indi¬ 
cates inverse. The prompted input in the 
following example is from the command 
completion user interface that I describe in 
the next section. 
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Enter distinct generators as a string 


e.g. RS means two generators R and S 
Generators: xy 

— Pressing 

ESC will abort the 

search 

Do you want these to generate the entire 

2 

group 

order 

= 

2 

X 

= 

B 

Y 

= 

B 

group? (y or n) Y 

5 

group 

order 

= 

4 

X 

= 

C 

Y 

= 

B 


8 

group 

order 

= 

6 

X 

= 

D 

Y 

= 

B 

Enter the exact order for each generator. 

13 

group 

order 

= 

8 

X 

= 

E 

Y 

= 

B 

Press Enter for no order specified 

18 

group 

order 

= 

10 

X 

= 

F 

Y 

= 

B 

X is of order 2 

22 

group 

order 

= 

12 

X 

= 

G 

Y 

= 

B 

Y is of order 

27 

group 

order 

= 

14 

X 

= 

H 

Y 

= 

B 


40 

group 

order 

= 

16 

X 

= 

I 

Y 

= 

B 

A relation is of the form LHS = RHS 

47 

group 

order 

= 

18 

X 

= 

J 

Y 

= 

B 

Put in LHS RHS or LHS ( if RHS is e ) 

52 

group 

order 

= 

20 

X 

= 

K 

Y 

= 

B 

<Press ENTER to quit> 

58 

group 

order 

= 

22 

X 

= 

L 

Y 

= 

B 


69 

group 

order 

= 

24 

X 

= 

M 

Y 

= 

B 

LHS RHS » xy y'x 

78 

group 

order 

= 

26 

X 

= 

N 

Y 

= 

B 

Generators: 

86 

group 

order 

= 

28 

X 

= 

0 

Y 

= 

B 

XY 

92 

group 

order 

= 

30 

X 

= 

P 

Y 

= 

B 

Orders: 

X= 2 

111 

group 

order 

= 

32 

X 

= 

B 

Y 

= 

D 


RELATIONS: 
XY 


Y' X 
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The group of quaternionic units (order 8) 
can be presented as a group with three 
generators, a, b, c, and relations ab = c, 
be = a, ca = b. The following chart shows 


the timing for the SEARCH command using 
different Forth implementations and com¬ 
puters. 


Timing for Search 

Groups32 ver 7.0 // Win32Forth high level 

1308270 ms 

Groups32 ver 7.0 // Win32Forth assembly language 

51300 ms 

Groups32 ver 7.0 // SwiftForth 

307656 ms 

Groups32 ver 6.4g // Gforth LINUX 

1794000 ms 

Groups32 ver 6.4g // Gforth on Internet 

2700000 ms 


The entry marked Win32Forth assembly 
language replaces just nine words by as¬ 
sembly language equivalents - the in¬ 
crease in speed is remarkable. 


The entries for Gforth compare a LINUX 
version of Forth on a personal computer 
with the timing for the same version (run¬ 
ning on a SUN server) accessed via Telnet 
on the Internet. 
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Notice that the fastest timing (about 1 
min) is for a personal computer using as¬ 
sembly language coding of a few time- 
critical commands, and the slowest, using 
the Internet, is about 45 min. With the ex¬ 
ception of the entry marked "Internet", the 
computer used has a Pentium III processor 
running at 450 Mhz. 

Return to Table of Contents 


The process of searching through a collec¬ 
tion of groups to find those with a given 
set of generators and relations is inher¬ 
ently time-consuming. This does not pre¬ 
clude the possibility that a better algorithm 
will make a major improvement. 
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User Interface 


When I first circulated Groupsl6 for review 
as an illustration of programming method¬ 
ology, several reviewers suggested that it 
might make a good instructional program 
- provided it could be made easier for stu¬ 
dents to use. 

Groups32, as presented so far, requires 
the user to learn a collection of command 
names, together with the parameters (and 
order) required for their use. As with simi¬ 
lar research systems, there are actually 
only a handful of top-level commands in 
the system. (The others are useful mainly 
for those who wish to write programs to 
extend the system.) It is not unrealistic to 


use the system with a "command line" in¬ 
terface, but there is a learning curve in¬ 
volved, making casual use more difficult. 
Typing commands and their arguments is 
slow. Also, a more constrained type of in¬ 
put would decrease the need for error 
trapping. 

I explored several methods to make these 
top-level commands more accessible and 
to ensure valid input. In one attempt, 
evoking a top level command provided 
prompts for parameters. In another, I pre¬ 
sented the top-level commands in a menu 
- clicking on a command initiated a 
prompt for input. 
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The approach I finally used, a "command 
completion" interface, appears to be the 
quickest and easiest to use. In this inter¬ 
face, the user is presented with a list of 
available commands: 


CENTER 

CENTRALIZER 

CHART 

CONJ-CLS 

COSETS 

EVALUATE 

GENERATE 

GROUP 

HELP 

ISOMORPHISM 

LEFT 

MENU 

NORMALIZER 

ORDERS 

PERMGRPS 

POWERS 

QUIT 

RESULT 

RIGHT 

SEARCH 

STOP 

SUBGROUPS 

TABLE 



The user types enough letters of a com¬ 
mand to distinguish it from the others. To 
generate a table, for example, the user 
needs only to type "T", since this is the 
only command that starts with T. To obtain 
a list of conjugacy classes, the user types 
"CON". The interface discards all letters 


typed beyond the bare minimum (so the 
user need not worry about mistyping sub¬ 
sequent characters). Once the interface 
identifies a command, it prints the rest of 
the command and prompts for any argu¬ 
ments needed to carry out its operation. 

Gl» CONJ-CLS Group Number 8 
{A} {BC} {DEF} 

Typing "CON" distinguishes the CONJ-CLS 
command and asks for a group number. 
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Here is a sketch of how I implemented the 
interface: 

The commands and information about 
them are put into a binary tree. Each 
node in the tree has a name field (for a 
string of a maximum size), the address of 
the node to the left, the address of the 
node to the right, and the execution ad¬ 
dress of the word used to carry out the 
command. 

32 CONSTANT Max#Cmds 
24 CONSTANT CmdSize 
0 VALUE #Cmds 

CmdSize CELL + CONSTANT EntrySize 
MaxtCmds CONSTANT #Nodes 


CREATE 'Treel #Nodes NodeSZ * ALLOT ALIGN 

: Left ( n_addr — l_addr ) SZ + @ ; 

: Right ( n_addr — r_addr ) 

SZ + 1 CELLS + @ ; 

: Exec ( n_addr — ) 

SZ + 2 CELLS + @ EXECUTE ; 

: Left! ( x n_addr — ) SZ + ! ; 

: Right! ( x n_addr — ) 

SZ + 1 CELLS + ! ; 

: Exec! ( e_addr n_addr — ) 

SZ + 2 CELLS + ! ; 

: Name! ( $ n_addr — ) DUP SZ BLANK $!! ; 

: Leaf? ( n_addr — flag ) 

DUP Right 0= SWAP Left 0= AND ; 
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Each new command is added to the tree to 
the left of an existing node if the name 
comes before in lexicographic order, and 
to the right if the command comes after. 
The commands are introduced by the word 
>CMD, which is used in the form 

>CMD <menu_name> <word_to_execute> 

The word_to_execute contains prompts 
and input routines to get any parameters 
needed from the keyboard. There is also a 
provision for including descriptive "help" 
information that will be printed if the user 
types HELP before typing the name on the 
menu. 


Here is an example of how the Table 
command is made into a word %Table that 
prompts for input: 

: %Table 

Help: 

This prints a table for the group requested 
(and makes that the current group). Elements 
are represented by letters A to Z and the 
symbols [ \ ] A _ and 
Help; 

." for " ['] Get-Grp CATCH 0= 

IF CR Table THEN ; 

>CMD TABLE %Table 
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I have written a collection of input com¬ 
mands. For example, Get-Grp is the com¬ 
mand to get a group number. In all cases, 
the system waits for the user to type a 
string, which is put into the terminal input 
buffer (TIB) as if it had been typed at the 
keyboard in an interactive session. The 
input is then processed by Forth words 
that read the input stream. Also, the 
cursor position is saved at the start of in¬ 
put, and the cursor is returned to this 
position if there is an error in the input. 


Here are other input commands that can 
be used for writing the prompted input 
sections. 


2VARIABLE Save-Pos 

: Get-TIB ( — ) AT? QUERY AT 

>IN @ 0 WORD 
COUNT TYPE >IN ! ; 

: Get-Num ( — n ) ." Number " AT? 

Save-Pos 2! 

BEGIN Save-Pos 2@ AT 

Get-TIB BL WORD 
COUNT ?DUP 0= THROW 
NUMBER? IF DROP TRUE 

ELSE 2DROP BEEP FALSE 
THEN 

UNTIL ; 

: Get-Grp ( — n ) ." Group " 

BEGIN ['] Get-Num CATCH 
IF TRUE THROW 
ELSE G-OK 

DUP 0= IF CR 10 SPACES THEN 

THEN 
UNTIL ; 
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: Get-Set ( — set ) 

." set: { " Get-TIB ." }" 

0 0 WORD COUNT 2DUP UPPER 0 

?D0 DUP C@ ID - 

DUP 0 MaxOrd 1- BETWEEN 
IF ROT Member! SWAP 
ELSE DROP THEN 
1 + 

LOOP DROP ; 

: Get-Ele ( — ele ) ." element " 

0 BEGIN 

DROP KEY UPC ID - 
DUP 0 Gord 1- BETWEEN 
DUP 0= IF BEEP THEN 
UNTIL ; 

: Get-Ele. ( — ele ) Get-Ele DUP .Ele ; 

Here is how the interface works. The char¬ 
acters typed by the user are added to an 
incomplete command string. Each time a 


character is added, the tree is searched to 
see if the string matches the leading part 
of a name in the tree. If there is no 
match, the user is informed (by a beep), 
and the latest character is not accepted. If 
there are multiple matches, the system 
waits for further characters. If there is only 
one match, the remainder of the name is 
typed and the stored command is exe¬ 
cuted. 

Since the user must type additional infor¬ 
mation from the keyboard, this interface 
has proved to be faster than one in which 
the user selects commands with a mouse. 
An additional feature of the interface is 
that it is very easy to add additional com¬ 
mands - so the interface is extensible. 
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It is easy to learn to use the interface, and 
no programming is required. As a result, 
Groups32 can be easily integrated into 
courses in abstract algebra. 

The original version of the interface used a 
single ordered list of commands. In later 
revisions I used a single binary tree. When 
I added permutations, I revised the code 
for the interface to enable using more than 
one tree. In Groups32 the permutation 
commands are in a second tree and there 
is a separate menu: 


CREATE ELEMENTS HELP INSTALL 

MAIN MENU MULTIPLY QUIT 

The PERMGRPS command in the main 
menu switches to this package, and the 
command MAIN here switches back to the 
main system. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 164 


Maintaining and extending the system 

Since a Forth application is extensible and 
Forth compiles incrementally, it is neither 
necessary nor desirable to make modifica¬ 
tions or extensions by revising and recom- 
piling the entire system. _ 


I emphasize this point because it is differ¬ 
ent from the situation that arises with con¬ 
ventional compiled languages. 

It has been found best, once a basic sys¬ 
tem is in place, to treat it as the imple¬ 
mentation of a language: The language 
should be documented and not frequently 
subjected to total revision. Additions and 
modifications can be made as extensions 
to the existing system. 


Most Forth implementations, including 
Win32Forth, have a provision to "save the 
system". The result is an executable that 
combines the underlying Forth system with 
the extensions added by the user. The 
combined system is then used as if Forth 
came with the user's words. 

When I have used Groups32 for instruc¬ 
tion, I have installed the executable on a 
file server, and the basic system is un¬ 
changed once installed. It is nevertheless 
possible to make changes and additions as 
the course progresses: Win32Forth reads a 
configuration file after the system is 
loaded. (Groups32 is saved as an executa¬ 
ble with the configuration file designated 
Groups32.cfg.) 
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The line I NCLUDE <pathname> updates, f is 

placed in the configuration file. The file 
"updates.f" is placed in a public directory, 
and the instructor has read/write privi¬ 
leges. All extensions, modifications, bug 
fixes, etc. are made in the file updates.f. 

Extensions 

Additions to the system can be made eas¬ 
ily in updates.f. For some classes I have 
added, for example, an additional package 
with number theory commands. Additional 
commands can be added easily to the 
menu in the commands completion inter¬ 
face. First create a version of the new 
command that prompts for its arguments 
and contains help information, then install 
it on the menu by 


>CMD <menu_narTe><:om[Tand_namE>. 

If this is done in updates.f, the new com¬ 
mand will appear on the menu whenever 
the system is loaded. 

Redefining Words 

It is sometimes desirable to change the 
action of an existing word. If a word is re¬ 
defined (i.e., a new definition is made us¬ 
ing the same name as a previously defined 
word), the new word is added to the end 
of the dictionary, and a dictionary search 
will find it first. Subsequent use of the 
name will find the new definition, but any 
word compiled with the older definition will 
continue to use the older action. 
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If >CMD is used with the new action, the 
command completion interface will use the 
new action. 

Most Forth systems provide a warning if a 
word is re-defined. 

: myword ; 

: myword this one prints" ; 

MYWORD isn’t unique 

The message "isn't unique" is not an error 
message - it's just a reminder to the user 
that a name is being used again. If this 
was done unintentionally, it can be undone 
in many Forth systems with FORGET (us¬ 
age: FORGET <wordname>), a word that 
will remove <wordname> and all subse¬ 
quent words from the dictionary. 


A new word can reference an older version 
in its definition. 

: print-one 1 . ; 

: print-one Here is one: " print-one ; 
PRINT-ONE isn't unique 
print-one Here is one: 1 

Forth does have a mechanism for recur¬ 
sion, but self-reference is not it. The cur¬ 
rent word is deliberately excluded from the 
dictionary search. In the example above, 
the new print-one uses the earlier version 
in its definition. 
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Patching 

As noted already, redefining a word does 
not change the existing action for words 
already compiled. It is sometimes neces¬ 
sary to change the behavior of a word in 
such a way that any words already com¬ 
piled use the new action. This can be done 
in at least two ways. 

1. Vectored Execution 

If it is known that a word may need to 
have its action changed, there is a mecha¬ 
nism for "deferred execution". The method 
defines a word that executes a stored ad¬ 
dress (vector). Changing the stored ad¬ 
dress changes the behavior of the word. 


I used this feature in the permutation 
package to permit a switch, after the fact, 
between left-to-right and right-to-left mul¬ 
tiplication: 

DEFER Direction 

: Left->Right [ 1 ] NOOP IS Direction ; 

: Right->Left [ 1 ] SWAP IS Direction ; 

The behavior of the word Direction is 
switched between a NOOP (no operation) 
and SWAP. This switchable word is used 
at the start of the code for P*. 
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Vectored execution is also used in the 
"print to file" code. Win32Forth, like many 
other Forth systems, uses deferred execu¬ 
tion for basic output words EMIT, TYPE, 
and CR. We can switch the behavior of 
these words. During printing to a file we 
install versions that not only print to the 
console but also send output to a file. 

Vectored execution can also be used effec¬ 
tively in development to test several ver¬ 
sions of a word: The code is written with a 
vectored word, and alternatives can be in¬ 
stalled "on the fly". 


For clarity of code, it is important to 
document, at the point where a vectored 
word is created, the type of behavior that 
will be filled in. 

2. Patching a definition 

The dictionary entries themselves occupy 
accessible memory. One can often change 
the behavior of an existing word (e.g., to 
fix a bug) by altering an existing dictionary 
entry directly. Details of this process de¬ 
pend on the implementation of Forth and 
are beyond the scope of this article. 


Return to Table of Contents 
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Summary 


The methodology I have described in this 
article involves creating a special-purpose, 
problem-oriented language. I use a base 
language (Forth) that permits users to get 
involved in the creation of the program¬ 
ming language. Forth provides access to 
memory, the input stream(s), and the 
process of compilation. 12 With this access 
one can implement almost any type of 
data structure or language syntax. 13 


12 Some Forth implementations also provide direct ac¬ 
cess to hardware - a feature used in instrumentation and 
control applications, but not in this article. 

13 Most of the language development in the case of 
Groups32 has been in the representation of data, the in- 


A special-purpose language provides a 
means of expression that is conceptually 
closer to the application area. In Groups32 
I have used the ability to conceptualize in 
terms of groups, subgroups, sets, etc. No¬ 
tice that, while early definitions are in 
terms of integers and addresses, later 
definitions are expressed more in terms of 
group elements, subgroups, etc. 


troduction of new data types, and the introduction of ap¬ 
propriate operations. 
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I have used a "middle out" approach to the 
design and implementation of a software 
system: A rudimentary, functional, system 
is created first. This core system is the 
major tool in the development of more ex¬ 
tensive and refined systems. 

I have used this methodology to produce 
systems for research as well as for instruc- 

Return to Table of Contents 


tion. These systems provide an interactive 
computing environment. They are exten¬ 
sible: In research, a system can be ex¬ 
tended and modified to meet the needs 
that arise as a project progresses; in in¬ 
struction, a system can be tailored to the 
needs of a particular class. 
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Resources 


1. The Forth implementation used in this 
work is Win32Forth, written by Tom 
Zimmer and Andrew McKewan. It is in 
the public domain. A complete ver¬ 
sion of Win32Forth and supplemen¬ 
tary materials can be found on the 
Forth Interest Group Web site, 
http://www.forth.org . 

2. A variety of non-commercial Forth 
implementations and other Forth re¬ 
sources can be found on the Forth In- 


Forth Interest Group Web site. There 
are also links at this site to the web 
pages of vendors of commercial Forth 
systems. 

3. More Forth resources are listed on my 
Web site. 

The following papers give background 
information on the approach to system 
development used in this paper, and 
source code for some of the features. 
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Supplementary Information 


This section contains the targets of hyperlinks. 
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The Forth Language 


Forth was invented (or perhaps discov¬ 
ered) during the 1970s by Charles Moore. 
Its original use was for applications involv¬ 
ing hardware control and instrumentation. 
It was used as a high-level, portable re¬ 
placement for assembly language pro¬ 
gramming. The language was designed for 
simplicity, flexibility, and speed. It is cur¬ 
rently used to assist in the design and 
testing of experimental hardware control 
applications. 

A major use of Forth is to program "em¬ 
bedded systems", i.e., computers that are 
included in other products. 


In many respects the advantages of Forth 
for developing experimental hardware con¬ 
trol systems are similar to the advantages 
we exploit for work in mathematics. It 
provides a simple comprehensible founda¬ 
tion to which a variety of features can be 
added. It is a language suited for building 
language features. 

Here are some of the attributes that make 
Forth a good choice for this work: 
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1. Forth is conceptually simple. It is fairly easy 
for a non-specialist to learn. 

2. Forth is suited for creating special-purpose 
languages. 

3. Forth provides an interactive computing 
environment. 14 

a. This speeds the development and testing 
of code. 

b. Programs can be easily written for interac¬ 
tive use. 

4. Forth programming creates a dictionary of 
useful words. This fine-grained modularity 
improves clarity and contributes to produc¬ 
ing reusable code. 


5. It is easy to introduce new types of alge¬ 
braic objects and operations seamlessly. 

6. Forth is small and efficient. 

7. Programs written in Forth execute at a 
speed comparable to conventional com¬ 
piled languages. 

8. Forth is extensible—extensibility includes 
introducing new data types and new lan¬ 
guage features. 

9. Details of the language implementation are 
accessible to the user. 

10. Forth supports many programming styles 
and paradigms 


14 By this we mean that a user's work with a system is a "session": commands can be issued from the keyboard to act on 
persistent data. The user can issue a command, see the result, and continue operations on the data. The systems 
produced by the methodology in this article also include the ability to create new commands during a session. 
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IF and THEN 


This level of detail is beyond the scope of 
my intended simple introduction. However, 
you may be curious about exactly how one 
can define words that do what we have as¬ 
serted that IF and THEN do. 


Here are the simplest forms of 
tions of IF and THEN used in 
Forth implementations: 

: IF COMPILE 7BRANCH HERE 0 , ; 

: THEN HERE SWAP ! ; IMMEDIATE 


1 F 

Start definition of IF 

COIVRI LE ? BFtANCH 

Compile conditional branch instruction 

HERE 

Put current address of dictionary pointer on stack 

O , 

Compile a zero as a placeholder to be filled in by 
proper address. 

' 

End the definition 

1 MVEDI ate 

Tag this as an immediate word 


the defini- 
traditional 

IMMEDIATE 
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THEN 

Start definition of THEN. Remember that IF 
left an address on the stack. 

HERE 

We put where we are now on the stack 

SVNAP 

Put location of missing address on top of current 
address 

! 

Store current address. 

; 

End the definition 

1 MVEDI ate 

Tag this as an immediate word 


In most Forth implementations additional information is put 
on the stack to check that conditionals are properly paired 
(an IF with a THEN, etc.) 
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Transporting Information 


Here, at the very lowest level, are exam¬ 
ples of the usefulness of the system and of 
the ability to program it. 

I have used a modified version of the TA¬ 
BLE word to export the group tables (to¬ 
gether with additional information) in a 
form suited for use by a Java Applet (see 
Groupsl5 J. 

I have mentioned the use of Groups32 to 
print information for input into spread¬ 
sheets for sorting an analysis. This was 
used by Evelyn Manalo (see Honors The¬ 
ses ) to develop a Maple procedure for 
printing the cycle decomposition of an abe¬ 
lian group and for work on the isomor¬ 
phism problem for groups of order 1-32 


(i.e., to find an efficient way to classify 
groups of order 1-32). 

Recently Groups32 was used to generate 
data suitable for use in Maple for a student 
research project on Cayley graphs (see 
Honors Theses ). The student, Sonja Willis, 
wanted to use the Maple networks package 
to print Cayley graphs with colored arrows 
rather than the green lines that the pack¬ 
age normally produces. To do this, one 
needs group data with the edges produced 
by the various generators kept separate. 
Some programming in Maple was needed 
to produce the colored arrows. The data 
were produced by Groups32. 







JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 179 


Here is an example for group #8 (S 3 ). The 
information transported was produced by a 
word PrintMaple (see below). 


8 Subgroups 

* = Normal subgroup 
Generators Subgroup 


0 

{ 

} 

*{ 

A 

} 


1 

{ 

D } 

{ 

A 

D 

i 

2 

{ 

E } 

{ 

A 

E 

i 

3 

{ 

F } 

{ 

A 

F 

i 

4 

{ 

B } 

*{ 

A 

B 

c 

5 

{ 

B D } 

*{ 

A 

B 

c 



J 5 


In this version, the generators were taken 
from the Subgroups command. B(l) and 
D(3) are generators. The Maple command 
generated by PrintMaple is a list consisting 
of the set of vertices, the number of gen¬ 
erators, and a list of sets of edges. 
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8 PrintMaple 

This produces the following output: 

Cay[8] := [{$0. .5},2, [\ 

{ [0,1], [1,2], [2,0], [3,4], [4,5], [5,3]\ 
}\ 

,{ [0,3], [1,5], [2,4], [3,0], [4,2], [5,1]\ 
}\ 

] ] : 

The key to understanding this is to realize that 
we have created a Forth word called PrintMaple 
that takes a group number as input and whose 
action is to print a Maple command. The de¬ 
tails follow on the next three pages. 
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: Out$ ( addr cnt — ) 

0 ?D0 

DUP C@ DUP 

ASCII # = 

IF DROP SWAP 0 .R 

ELSE EMIT THEN 

1 + 

LOOP DROP ; 

Making output strings with inserted numbers. 

Here is a version taking the address and count of a 
string as input. When a # occurs in the string, it prints 
the number on top of the stack. 

?\ GETXY DROP 40 > 

IF . " \ " CR 8 SPACES THEN ; 

This prints a Maple continuation character and line 
break if the line gets too long. 

VARIABLE PMFirst 

VARIABLE GenFirst 

These variables are used to determine whether an 
item is the first one being printed. 
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: PrintMaple ( grp — ) 

This prints the function header and set of vertices. 

>Group CR 

The first # is the number of the group and the sec¬ 

Gord 1- Gnum 

s" Cay [#] := [{$0. .#}, " 

Out$ PMFirst ON 

ond is Gord-1 

Gnum 144 < 

Generators are obtained from the Subgroups 

IF Gnum Generate-Subgroups 

command. The last subgroup generated is G itself. 

#Subgroups 1- Generators 

Groups32 does not compute subgroups for group 

ELSE 2198 \ {bcehl} 

THEN 

number 144 - so this is treated as a special case. 

DUP #Set S" #,[\" Out$ 

CR 8 SPACES 

Print number of generators 





JOMA 


Journal of Online Mathematics and its Applications 


Wavrik- 183 


GenFirst ON 

32 0 DO I OVER MEMBER? 

IF \ Output for a generator 
GenFirst @ 

IF {" GenFirst OFF 

ELSE ." ,{" THEN 

Test if loop index is a generator. 

We will print a set bracket { only before the first set 
of edges, otherwise we print,{. 

PMFirst ON 

For-All-Elements 

DO 

PMFirst @ 

IF [" PMFirst OFF 

ELSE , [" THEN 

For each generator we print the corresponding 
edges [a,ag]. We print a comma between edges. 

I J G* I S" #, #] " Out $ ?\ 
LOOP 

}\" CR 8 SPACES 

THEN 

LOOP DROP . " ]] : " 

CR ; 

We print the pair for an edge and a continuation 
character if necessary. Then close the set of edges 
with}. 

Finally, end the structure. 
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Tables for Abelian Groups 

The permutation package can be used to 
produce a group table for any abelian 
group. The Fundamental Theorem of Abe¬ 
lian Groups says that any finite abelian 
group is a product of cyclic groups. The 
permutation package allows us to create 
products of cyclic groups by entering dis¬ 
joint cyclic permutations for generators. 
Thus Z 2 x z 4 can be realized as the sub¬ 
group of S 6 generated by ( 1 2 3 4 ) and 
(5 6). This group is produced in the per¬ 
mutation package by: 

6 Size! 

Gen: (1234) 

Gen: ( 5 6 ) 

Make-Group 

Show-Table 


1 2 

2 1 

3 4 

4 3 

5 6 

6 5 

7 8 

8 7 


3 4 

4 3 

5 6 

6 5 

7 8 

8 7 

1 2 

2 1 


5 6 

6 5 

7 8 

8 7 

1 2 

2 1 

3 4 

4 3 


7 8 

8 7 

1 2 

2 1 

3 4 

4 3 

5 6 

6 5 


To generate the tables for abelian groups 
of orders 17-31, I automated this idea. I 
actually generated the tables with an ear¬ 
lier version of the permutation package, 
which used a slightly different representa¬ 
tion for the permutations and which di¬ 
rectly stored the groups in the internal 
representation. Here is a more modern ad¬ 
aptation that introduces an interesting 
useful possibility: Generate a script file, 
and then run the script file to produce the 
desired results. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik- 185 


Script Files 

A script file is a file of commands. It 
could be a file of Forth/Groups32 com¬ 
mands (as in this case). It could also be a 
file of commands for some other system - 
we often produce script files for Maple. 
The easiest way to make a script file is to 
capture output to the screen in a file. I 
wrote for Groups32 a "print to file" utility 
for this purpose. A script file is loaded into 
Forth by INCLUDE <filename>. The Forth 
system then reads and executes com¬ 
mands from this file as if they were being 
typed at the keyboard. 

First we need a Groups32 command called 
CycDecomp. The stack diagram is ( si s2 
... sk n — ), where si, s2, ... , are the sizes 
of the component cycles, and n is the sum 


of the sizes. (It will determine the Sn in 
which they are considered to be permuta¬ 
tions.) The output of CycDecomp will be 
the commands needed by the permutation 
package to produce the table of the asso¬ 
ciated group. 

235 CycDecomp 

5 Size! 

Gen: (123) 

Gen: ( 4 5 ) 

Make-Group 

Show-Table 

Notice that the word CycDecomp does not 
execute these commands - it just prints 
them on the screen. If the print-to-file 
feature is active, these commands are also 
captured to a file. 
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Win32Forth can also print the console win¬ 
dow to a file - so, for a limited amount of 
output, it may be sufficient to use this fea¬ 
ture. 

The abelian groups of order n correspond 
to partitions of n. Thus, for example, we 
can partition 6 = 3 + 2 + 1. If we use 
1 2 3 6 CycDecomp, we obtain commands 

6 Size! 

Gen: (123) 

Gen: ( 4 5 ) 

Gen: ( 6 ) 

Make-Group 

Show-Table 


Executing these commands produces the 
table 

1 2 3 4 5 6 

2 1 4 3 6 5 

3 4 5 6 1 2 

4 3 6 5 2 1 

5 6 1 2 3 4 

6 5 2 1 4 3 

By inspection, this table represents an 
abelian group of order 6. "1" is the iden¬ 
tity. Element "4" has order 6 - so this is 
the table of Z 6 . We must, therefore, select 
from all partitions of n those that produce 
non-isomorphic groups. Here are the 
commands for groups of orders 17-22: 
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17 

17 

CycDecomp 

2 

4 

3 

9 

CycDecomp 



29 

29 

CycDecomp 

2 

3 

3 

8 

CycDecomp 


8 

3 

11 

CycDecomp 

2 

3 

5 

10 

CycDecomp 


2 

9 

11 

CycDecomp 


5 

5 

10 

CycDecomp 



31 

31 

CycDecomp 



19 

19 

CycDecomp 



25 

25 

CycDecomp 

2 2 2 

2 

2 

10 

CycDecomp 

2 

2 

5 

9 

CycDecomp 


2 

13 

15 

CycDecomp 

2 2 

2 

4 

10 

CycDecomp 


4 

5 

9 

CycDecomp 

3 

3 

3 

9 

CycDecomp 

2 

4 

4 

10 

CycDecomp 


3 

7 

10 

CycDecomp 


3 

9 

12 

CycDecomp 

2 

2 

8 

12 

CycDecomp 


2 

11 

13 

CycDecomp 



27 

27 

CycDecomp 


4 

8 

12 

CycDecomp 



23 

23 

CycDecomp 

2 

2 

7 

11 

CycDecomp 


2 

16 

18 

CycDecomp 

2 2 

2 

3 

9 

CycDecomp 


4 

7 

11 

CycDecomp 



32 

32 

CycDecomp 


The definition of Show-Table was changed, 
for the purpose, to print the tables in the 
format used to load tables from text files. 


Here is the format: 


7 ->Group 

0 7 >ROW 0123456 

1 7 >ROW 1234560 

2 7 >ROW 2345601 

3 7 >ROW 3456012 

4 7 >ROW 4560123 

5 7 >ROW 5601234 

6 7 >ROW 6012345 


Here is the procedure: 

1. Produce a word CycDecomp to make en¬ 
tries in a script file containing commands 
for the permutation package. 

2. INCLUDE the script file and produce tables 
for the groups in a desired format. This is 
saved as a second script file of Groups32 
commands for loading groups. 

3. INCLUDE the second script file to load the 
groups. 
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G* and Assembly Language 

Here are the three definitions for G* that 
were used in the timing chart: 

1. Groups 16 (1990) where elements are represented by A, B, 
C„.. .(but BLOCK has been redefined to use memory) 

: 'ELE ( elel ele2 - addr ) 

\ compute address of product 

ID - SWAP ID - GORD * + 

\ offset from base address 

GBLK BLOCK GOFFS + 

\ address of start of table 

+ ; \ address of element 

: G* ( elel ele2 - elel*ele2 ) 

'ELE C@ ; 


2. Groups32 - elements are represented by 0, 1,2,... 

: G* ( elel ele2 - product ) 

MaxOrd * + \ compute offset of product 

Grp @ + \ add the base of the table 

C@ ; 

3. Groups32 with all tables 32 x 32, elements repre¬ 
sented by 0, 1,2,..., and code written in assembly 
language for the Intel 80x86. 

CODE G* ( i j — i* j ) \ ASSUME MaxOrd=32 

SHL EBX, 5 \ Multiply j by 32 

POP EAX ADD EBX, EAX \ Add i 
MOV EAX, Grp [EDI] \ Add group offset 
ADD EBX, EAX 

SUB EAX, EAX \ byte at this address 
MOV AL, [EBX] [EDI] 

MOV EBX, EAX 
NEXT END-CODE 

There are 9 words in Groups32 that are 
coded in assembly language. 
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WORD 

Most Forth words take parameters from 
the stack and return results to the stack. 
The Forth word WORD reads the input 
stream. The object is to convert the text 
following it into an internal form (a 
counted string). 

A counted string is stored in memory as a 
count (one byte) followed by the ASCII 
codes of its characters. The string ABC 
would be stored as 


3 

65 

66 

67 


Once it is stored, the address of the count 
byte is passed to other words on the stack. 
We imagine that a string is on the stack. 
We say to ourselves "a string is on the 
stack". We tell our friends "a string is on 
the stack". In reality, the stack contains 
just the address of the count byte, but we 
can conceptualize as if the stack contained 
the string itself. This is typical of the way 
that we make our system appear to handle 
sophisticated types of data on the stack. 
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WORD ( delim - c-addr ) is an input word. 
It reads the text that follows it when it 
acts. It stores the text as a counted string 
in a temporary location (in most Forth im¬ 
plementations, this is at the end of the dic¬ 
tionary), and it returns the address of the 
counted string. WORD is usually used in 
definitions that immediately do something 
with the string (process it or move it else¬ 
where). 

WORD is the word the system uses when it 
interprets the input stream word by word. 
It could, but usually doesn't, use a tempo¬ 
rary storage mechanism. Thus each in¬ 
coming word is placed in the same tempo¬ 
rary location. Forth gets away with this be¬ 
cause the string is immediately processed. 


If you try to see what WORD does (as a 
referee of this paper did) by typing a 
command line, you will find it is a shy little 
animal: Just trying to look at it chases it 
away. You might try 

BL WORD ABC COUNT TYPE 

and you will get TYPE as the response. 
What happens is that BL WORD does read 
the ABC and make it a counted string 
(stored at the end of the dictionary). But 
then the Forth interpreter uses WORD 
again to read the string COUNT (and 
stores it at the end of the dictionary). It 
looks this up in the dictionary and exe¬ 
cutes it. It then reads TYPE (and stores it 
at the end of the dictionary) and looks it 
up and executes it. 
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The ABC string was overwritten by COUNT 
and then by TYPE - all of which were put 
at the same address. The very act of try¬ 
ing to observe the result of BL WORD ABC 
destroys what you are trying to observe! 

This is really not as horrible as it seems. 
The counted string produced by WORD is 
put in a temporary place known to be 
available - so everybody uses that place. 
If you really want to read the input stream 
yourself, you can - but you must do some¬ 
thing with the counted string the system 
uses WORD to read the next command. 


My suggestion, if you want to make a 
string package, is that you make some 
sort of temporary storage mechanism. I 
personally use the method discussed in 
this article. Others have used a string 
stack or a circular buffer. Then make an 
input word that uses WORD but immedi¬ 
ately moves the string to a temporary 
storage location. 

It is also legal to use WORD and immedi¬ 
ately do something with it. 




ML* **##*#****« 
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Virtual Memory and >GROUP 


Early versions of Forth ran on machines 
using 16 bits for integers and addresses. 
This meant that one could address only 2 16 
= 65536 memory locations. This 64k of 
address space included the code for the 
operating system and the language im¬ 
plementation, leaving very little space for 
the application code and data. Virtual 
memory is a scheme for treating the disk 
(it was floppy at the time) as if it were real 
memory: Data could be stored on disk 
and read to "real" memory when needed - 
and unneeded data would be paged out 
from "real" memory to disk. The trick was 
to make this process of swapping between 
disk and real memory as transparent as 
possible. 


Forth was developed at a time when com¬ 
puters were (by present standards) very 
small. An entire University was served by a 
computer with 4 Meg of RAM. A personal 
computer (such as my first one) was con¬ 
sidered rather powerful with 16k of RAM. 
When I expanded to 64k I was asked, se¬ 
riously, what I would do with all that 
memory. 

The operating systems of the time were 
too small and unsophisticated to include a 
virtual memory system. Forth was re¬ 
garded as rather advanced because it had 
such a scheme. 



JOMA 


Journal of Online Mathematics and its Applications 


Wavrik - 193 


The virtual memory scheme works like 
this: Imagine the disk to be organized as 
"blocks" of 1024 bytes. We number the 
blocks: 

DISK: 0 12 3 4 


MEMORY: 

A(n)u? B(n)u? 


Real memory has room for only two boxes 
that can store data read from disk. There 
are two vital pieces of information stored 
about these two "buffers": the number of 
the disk block currently in memory, and a 
flag (true or false) to indicate whether this 
information has been "updated" or 
changed. 

If you want to use an item stored in a disk 
block, it must be read into one of the buff¬ 
ers. It can then be read and changed just 
as if it is in memory. If it has been 
changed, it should be saved to disk before 
the buffer is reused. All this is done by the 
word BLOCK. 
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BLOCK ( n - addr ) 

IF 

block n is already in a buffer, return its ad¬ 
dress 

ELSE 

get the least recently used buffer 

if it is marked "update", write its contents to 
disk 

read the data from disk block n to this 
buffer 

mark the buffer as not updated 

return the address of the start of the buffer 

As you can see from this, if a certain disk 
block contains the data for one or more 
groups, these data are read into a buffer - 
but then all computations on that group 
can be performed using the data in mem¬ 


ory. It is almost as if the group data were 
in memory all the time. We just make all 
computations refer to the address (re¬ 
turned to us) of the start of the buffer that 
holds the data. 

So we have stored each group table in a 
certain disk block (the number is GBLK). 
The start of the table is a certain distance 
or offset from the start of the block 
(GOFFS). The rows of the table are stored 
one after another, so we need to know the 
group order (GORD) to find a particular 
product. We will also need to know the 
number of the current group (GNUM). 
This vital information about all of the 
groups was stored in block 0. 
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So the scenario is this: You want to study 
group n. You go to block 0 (IDXBLK) 
where the triples GORD, GBLK, GOFFS are 
stored. Each of these (in the original 16- 
bit version) occupies 2 bytes. To find the 
triple for group n, you must add 6n to the 
address of the start of the block. So, if 
Base denotes the address of the start of 
block 0, then 

GORD is at address Base + 6n, 

GBLK is at address Base + 6n + 2, 

GOFFS is at address Base + 6n + 4. 

To get the address of the base of the 
group table we do GBLK BLOCK GOFFS +. 
GBLK BLOCK will load the disk block con¬ 
taining the table into memory and return 
the base address of the block. Adding 


GOFFS will give the address of the start of 
the table. 

Once the triple GORD, GBLK, GOFFS is 
read from block 0, we no longer need to 
have block 0 in memory. When we do 
group computations in the current group, 
the data for that group are written from 
disk to memory and stay in position. No 
further disk access is necessary to com¬ 
pute in this group. 
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: >GROUP ( n — ) \ set GORD GBLK and GOFFS 

( n ) DUP TO GNUM \ save group number 

( n ) 6 * 

( 6n ) IDXBLK BLOCK + \ address of triple 

( addr ) DUP @ TO GORD 

( addr ) 2+ DUP @ TO GBLK 

( addr+2 ) 2+ @ TO GOFFS ; 


The original data files used 16-bit integers, 
so the three elements of the triplet occu¬ 
pied 2 bytes each. When an original data 
file is used with a 32-bit version of Forth, 


the fetch operation @ must be replaced by 
W@ so that it gets a 16-bit rather than a 
32-bit value. 
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Location in Table 

In Groupsl6 the table for a group of order 
m was stored as m x m consecutive bytes 
to save memory. In Group32 all tables are 
stored as 32 x 32, regardless of the group 
order, to increase speed. 

The group table for Z4, of order 4, looks 
like this in two dimensions: 


A 

B 

C 

D 

B 

C 

D 

A 

C 

D 

A 

B 

D 

A 

B 

C 


In memory this is stored by rows: 

ABCDBCDACDABDABC 

If we number the rows and columns 0, 1, 
2, 3, then the entry in row r and column c 
occurs 4r + c slots from the start of the 
table. Thus the first row (r=0) is at offsets 
0, 1, 2, 3 from the start of the table. The 
next row is at 4, 5, 6, 7. 

The word 'ELE ( elel ele2 - addr ) calcu¬ 
lates the address of the product of elel 
and ele2. Remember that in Groupsl6 the 
elements are represented as ASCII charac¬ 
ters - so we must subtract ID (the code 
for "A") to get r and c. 
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:'ELE ( elel ele2 - addr) 


ID - 

SWAP ID - 
GORD * + 
GBLK BLOCK 
GOFFS + 

+ ; 


( elel c ) 

( c r ) 

( r*ord+c ) 

\ address of block for table 
\ add the offset to get 
\ base address of table 
\ get address of element 


: 'ELE ( elel ele2 - addr ) 

Compute row and column numbers 
compute displacement of element 
(r*ord+c) 

get base address of stored table 
get address of element ; 
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Definition of Eorder 

Here is a more detailed discussion of the 
code. 

: EORDER ( x - order ) 

We plan to take the original element and create its pow¬ 
ers xn in a loop. We will quit when we get a power that is 
the identity element. 

ID 

We put the identity on the stack as x'D. So the stack has 
(xx"T>) 

GLIMITS DO (X X A n ) 

The loop is taken over all elements of the group. I should 
actually have made this GORD 0 DO - it would then have 
been clearer that we are looping over powers. I just didn't 
think of it at the time. The ( x xOn ) above is a comment. 
It is helpful to know that each time we are at the start of 
the loop this is what is in the stack. 

OVER 

If we had (xxdn ) in the stack, we now have (xxdn x). 


G* 

Now we have (xx^n+1]). 

DUP ID = 

Test if it is the identity but DUP it first so that it is not re¬ 
moved from the stack. 

IF 

2DR0P 

This removes the x and x^n+1]. 

I ID - 1+ 

Due to the way the loop limits were set up, this gives the 
exponent. 

LEAVE 

This causes an exit from the loop. 

THEN 

Notice that if we did not exit from the loop we now have 
( x x^n+l] ) in the stack, and the loop repeats with the 
stack as indicated by the stack diagram after DO. 

LOOP ; 

♦ ** * * ** * * * * ** * * *** * ** * * *9* * »* ■ »♦* * ** «- 
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Number of Squares 

CREATE SQ 20 ALLOT 
: #SQUARES ( grp# — n ) >GROUP 
SQ 20 ERASE 
GLIMITS DO 

I I G* ID - SQ + 1 SWAP C! 
LOOP 

0 

20 0 DO SQ I + C@ + LOOP ; 

This code is short, but might benefit by 
some analysis. The point here is to run 
through the group computing the squares 
of the elements. Each time a square is 
computed, a record is made that this ele¬ 
ment is a square. The structure is this: 


: #Squares ( grp# — n ) 

>Group 

Initialize-SQ 
Mark-Squares 
Count-Marks ; 

An array-like structure, SQ, is used to re¬ 
cord the marks. There is no standard AR¬ 
RAY word in Forth, but one can easily be 
defined (and a few are in the current 
Groups32). Rather than use general pur¬ 
pose machinery, we have used a more di¬ 
rect approach: 
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SQ will be the name for a block of 20 
bytes. When the name SQ is used it will 
return the address of the start of this 
block. So this word sets all these locations 
to 0: 

: Initialize-SQ SQ 20 ERASE ; 

Now these memory locations are intended 
to be associated with group elements. 
Since (in the 1990 version) we repre¬ 
sented group elements by the ASCII codes 
for A, B, C, .... , we must subtract ID to 
find the offset of a given element, x ID - 
SQ + will give the address in the SQ list 
corresponding to x. 1 <addr> C! stores 
the byte 1 at the given address. So now 
we can understand the code for marking 
the squares: 


: Mark-Squares 
GLIMITS DO 

I I G* \ compute the square 
ID - SQ + \ find address of the entry 
1 SWAP C! \ store 1 in this location 
LOOP ; 

After running this code, the SQ array will 
have a 1 in any place that corresponds to 
a square and a 0 otherwise. We therefore 
need only to run through the array and 
tally the number of l's. 

: Count-Marks ( — cnt ) 

0 

20 0 DO ( total ) 

SQ I + C@ \ gets the 1 or 0 
+ \ add it to the total 

LOOP ; 
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ZIP Files 


You can download a ZIP file containing executables (binaries) for the two versions of the 
groups system discussed. In most cases the ZIP file can be downloaded just by clicking 
here. 


In some cases it might be necessary to obtain the ZIP file using anonymous FTP. The pro¬ 
cedure for doing this is: 

1. Connect to ftp math.ucsd.edu 

2. Use "anonymous" as the username 

3. Use your email address as password 

4. cd pub/jwavrik 

5. type "binary" to ensure a binary transfer 

6. type "get wavrik_executables.zip" 






